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Abstract. We study the problem of conformal deformation of Riemannian structure to constant 
scalar curvature with zero mean curvature on the boundary. We prove compactness for the full set 
of solutions when the boundary is umbilic and the dimension n < 24. The Weyl Vanishing Theorem 
is also established under these hypotheses, and we provide counter-examples to compactness when 
n > 25. Lastly, our methods point towards a vanishing theorem for the umbilicity tensor, which is 
anticipated to be fundamental for a study of the nonumbilic case. 



1. Introduction 

The Yamabe problem consists of finding a constant scalar curvature metric g which is pointwise 
conformal to a given metric g on an n-dimensional (n > 3) compact Riemannian manifold M without 
boundary. This is equivalent to producing a positive solution to the following semilinear elliptic 
equation 

71 + 2 

(1.1) L g u + Ku^ = 0, on M, 

where K is a constant, L g = A g — c(n)R g is the conformal Laplacian for g with scalar curvature R g , 
and c(n) = gT^Ern • If u > is a solution of (11.11) then the new metric g = u™- 5 g has scalar curvature 
c(n)~ l K. This problem was solved in the affirmative through the combined works of Yamabe |46| . 
Trudinger (33], Aubin [5] and Schoen [39] (see also [27] for a complete overview). 

From an analytic perspective the Yamabe problem has proven to be a rich source of interesting 
ideas. The complete solution of the problem was the first instance of a satisfactory existence theory 
for equations involving a critical exponent, where the standard techniques of the calculus of variations 
fail to apply. When the first eigenvalue of the conformal Laplacian is positive, which is equivalent 
to the case of positive Yamabe invariant Y{M) (see [UJ for a definition), solutions to (II. ip are not 
unique, and it is known that the set of solutions can be quite large ([371 [36]). 

Therefore it becomes natural to ask what can be said about the full set of solutions to (jl.ip when 
Y(M) > 0. While this set is noncompact in the C 2 topology when the underlying manifold is S n with 
the round metric (see |36j), when M is not conformally equivalent to the round sphere compactness 
was established in various cases, namely by Schoen [38] in the locally conformally flat case, Schoen 
and Zhang [33] in three dimensions, Druet [15] for n < 5, Marques [32j for n < 7, Li and Zhang 
[29\ 150] for n < 11. However in a surprising turn of events, counterexamples to compactness were 
found by Brendle [7J when n > 52 and subsequently by Brendle and Marque [9J for 25 < n < 51. 
Finally, Khuri, Marques and Schoen [25] proved that compactness does hold in all remaining cases, 
that is, for n < 24. See [TO] for a survey of various compactness and non-compactness results for the 
Yamabe equation. 
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An obvious extension of such problems is to consider manifolds with boundary. In this case one 
would like to conformally deform a given metric to one which has not only constant scalar curvature 
but constant mean curvature as well. This problem is equivalent to showing the existence of a positive 
solution to the boundary value problem 



where v g is the unit outer normal and n g is the mean curvature. If such a solution exists then 

4 

the metric g = u n ~ 2 g has scalar curvature c(n)~ K and the boundary has mean curvature c. This 
Yamabe problem on manifolds with boundary was initially investigated by Escobar [171 H8j . who 
solved the problem affirmatively in several cases. With contributions from several authors (see 

[E3I331I3HI23I21IHI31IS1II2]), most of the cases have now been solved. 

Notice that if K ^ and c ^ then both the equation and the boundary condition are nonlinear. 
In order to simplify the problem, it is customary to assume then that one of them is linear, that is, 
that either K or c is zero. Geometrically, this corresponds to deforming the manifold to one with 
either constant nonzero scalar curvature and zero mean curvature on the boundary {K ^ 0, c = 0) 
or zero scalar curvature and constant nonzero mean curvature on the boundary (K = 0, c / 0). In 
this paper we will focus on the first of these two cases. 

In analogy to the case of manifolds without boundary, where the round sphere provides the canon- 
ical example of noncompactness, when the manifold has boundary and is not conformally equivalent 
to the round hemisphere, the question of compactness of solutions arises. Compactness was proven 
by Han and Li [23] when the scalar curvature is negative (K < 0) and the mean curvature is zero 
(c = 0), and also when the scalar curvature is positive (K > 0) with no restriction on the mean cur- 
vature but with the extra hypotheses that the manifold is locally conformally flat and the boundary 
is umbilic; by Felli and Ahmedou |20j when the scalar curvature is zero (K = 0), the mean curvature 
positive (c > 0), the manifold is locally conformally flat and the boundary umbilic (see also|21j); and 
by Almaraz [2] when the scalar curvature is zero (K = 0), n > 7, and a generic condition on the 
trace-free part of the second fundamental form holds. 

It is natural to consider subcritical approximations to equation (jl.2p . where a priori estimates are 
readily available. Thus we define 



from now on that K > 0. Then our main result may be stated as follows. 

Theorem 1.1. (Compactness) Let (M n ,g) be a smooth compact Riemannian manifold of dimension 
3 < n < 24 with umbilic boundary, and which is not conformally equivalent to the standard hemisphere 
(S+,go). Then for any e > there exists a constant C > depending only on g and e such that 



for all u G U 1+£<p< n+2 <3? p , where < a < 1. 

This theorem is established by a fine analysis of blow-up behavior at boundary points; such a 
fine analysis was carried out for interior blow-up points in [25]. The entire problem is reduced to 
showing the positivity of a certain quadratic form on a finite dimensional vector space, which may 
be analyzed in a similar manner as is done in the appendix of [25]. Of course this theorem also 
relies on the Positive Mass Theorem of General Relativity, in its usual form. That is, although we 



(1.2) 





C < u < C and \\ u ||c 2 .«(M)^ C 
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are concerned with manifolds having boundary, we are still able to use the standard Positive Mass 
Theorem by employing a doubling procedure. 

Another key feature of our approach is to employ a version of conformal normal coordinates 
adapted to the boundary, which elucidates the dependence of various geometric quantities on the 
conformally invariant umbilicity tensor and Weyl tensor. This coordinate system can be thought 
of as a good compromise between traditional conformal normal coordinates [27J and the so-called 
conformal Fermi coordinates [33]. This is because although the latter has been shown to be a powerful 
tool to study the Yamabe problem on manifolds with boundary a critical part of the compactness 
result in [25] is the proof of the positivity of the quadratic form mentioned earlier. This proof makes 
substantial use of the the radial symmetry coming from normal coordinates and we would like to 
preserve as much as possible of that original argument. 

In general, it is expected that wherever blow-up occurs, these conformally invariant quantities will 
vanish to high order because, up to a conformal change, the geometry of the manifold resembles that 
of a sphere near the blow-up. As we are assuming that the boundary is umbilic here, we focus on 
the Weyl tensor. In this regard we prove 

Theorem 1.2. (Weyl vanishing) Let g be a smooth Riemannian metric defined in the unit half n-ball 
Bi, 6 < n < 24. Suppose that there is a sequence of positive solutions {ui} of 



Pi € (1, ^z§]> such that for any e > there exists a constant C(e) > such that sup B +^ B + < C(e) 

and limj^oo (sup B + Uj) = oo. Assume also that Bf n M™" 1 is umbilic. Then the Weyl tensor W g 
satisfies 



Remark 1.3. It may appear that since the boundary is umbilic, the proofs of theorems 11.11 and 11.21 
should follow directly from |25] by applying a reflection argument. However, the techniques employed 
in [25] require a higher degree of regularity than what is typically available from a simple reflection 
of the metric. 

In analogy to the case without boundary, one wonders if theorem 11.11 is false when n > 25. We 
have also been able to answer this question. 

Theorem 1.4. Assume that n > 25. Then there exists a smooth Riemannian metric g on the 
hemisphere S™ and a sequence of positive functions m € C°°(S™), such that: 

(a) g is not conformally flat (so in particular (S™,p) is not conformally equivalent to (S+,go), 
where go is the round metric), 

(b) dS™ is umbilic in the metric g, 

(c) for each i, Ui is a positive solution of the boundary value problem: 




= 



in B±, 

on Bf n W 



Wg\( X ) < C\X 



for some integer i > 




n + 2 



where K is a positive constant, 

(d) supgn Ui — >■ oo as i — > oo. 
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Together, theorems 11.11 and 11.41 give a complete answer to the question of compactness of solutions 
to the Yamabe problem on manifolds with umbilic boundary in the positive scalar curvature setting 
(Almaraz has proven an analogue to theorem 11.41 for scalar-flat manifolds [1]). 

The proof of theorem 11.41 relies heavily on |9]. In fact, with [JJ 19] at hand, the idea to proof 
theorem 11.41 is not complicated. Brendle and Marques' construction is a perturbation of the round 
sphere (S n ,go). Although their solutions are constructed on S n rather than S 1 ", they "almost" satisfy 
the boundary condition. We can therefore slightly modify Brendle and Marques' solutions in order 
to produce a blow-up sequence for the hemisphere. 

One obvious consequence of theorem 11.11 is to give an alternative proof of the solution to the 
Yamabe problem, allowing us to compute the total Leray-Schauder degree of all solutions to (|1.2|) 
(with c = 0), and to obtain more refined existence theorems. This is discussed at the end of the 
paper (see section [1~4"|) . 

As mentioned earlier, certain conformally invariant quantities are expected to vanish to high order 
at a blow-up point. In particular such behavior is expected for the umbilicity tensor when the 
boundary is not umbilic. In this regard, we expect the following. 

Conjecture 1.5. Let g be a smooth Riemannian metric defined in the unit half n-ball Bf, 4 < n < 
24. Suppose that there is a sequence of positive solutions {ui} of 

= 0, in B+, 

on~BfnW l -\ 

Pi £ (1) ^rfL such that for any e > there exists a constant C(e) > such that sup B +^ B + Uj < C(e) 
and Hindoo (sup B + Uj) = oo. Then the umbilicity tensor T g satisfies 

\T g \(x) <C\x\ m , 
for some integer m > r± ^- . Moreover, if n > 6 we also have 

\W g \(x) < C\xf, x G B+, 

for some integer t > 

Proving this conjecture would be a key step towards a compactness theorem for manifolds with 
non-umbilic boundary. In fact, one of the main ingredients of our proofs is to estimate several relevant 
quantities in terms of the umbilicity tensor and its derivatives at the origin. The vanishing of these 
terms should allow one, at least in principle, to adapt the ideas presented here to the non-umbilic 
case. 

2. Setting, notation, and basic definitions 

Let M n be a n-dimensional Riemannian manifold with smooth boundary, and let {gi}^Zi be a 
sequence of metrics on M converging in C k (M) to a metric g, where k is large and depends only on 
n. Let {iti} be a sequence of positive solutions of the boundary value problem 




(2.1) 



' L gi Ui + Kf7^ u f = 0, inM, 
B gi Ui = d Vg _Ui + ^K gi Ui = 0, on dM, 



where L gi = A gt — c(n)R gi , c(n) = 4 ^_ 2 -^ , R gi is the scalar curvature of the metric <?j, K = n{n — 2), 
v gi is the outer unit normal, K gi is the mean curvature of the boundary, {/j} is a sequence of smooth 
positive functions converging in C 2 (M) to a smooth positive function /, 1 < pi < Si = —p%. 
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Lg is referred to as the conformal Laplacian, and the boundary value problem (J2JJ) is conformally 
invariant (see proposition IA.2|) . 

In conformal normal coordinates (see [271 El E] ) centered at a point p, we write g(x) = exp(/i(x)), 
where h is a smooth function taking values in the space of symmetric nxn matrices. From standard 
properties of conformal normal coordinates it then follows that x^hij(x) = 0, and trhij(x) = 0(r N ), 
where r = dist 9 (x,p) and iV is arbitrarily large. We also have det gij = 1 + 0(r N ). 

In most of the text we will identify the center p of normal coordinates with the origin. We will 
write Ui(x) instead of Uj(exp p (x)) and \x\ instead of dist g (x,p). Since N in the above expressions is 
as large as we want, we will often ignore the 0(r N ) contribution in the volume element and write 
dvolg(x) = dx. 

The proofs of theorems 11.11 and 1 1 . 2 1 dep end crucially on finding a good approximation to the scalar 
curvature in terms of polynomials. To this end we define, in conformal normal coordinates 

(2.2) H ij (x)= 

2<H<n-4 

where hi^ a are the coefficients of the Taylor polynomials centered at the origin. Notice that we will 
sometimes use jQ to denote Taylor coefficients at the origin — which are multiples of derivatives 
evaluated at the origin rather than the derivatives themselves. 

We then have = + 0(|x| n_3 ), = Hji, x : 'H i j(x) = 0, and ti Hij(x) = 0. Put also 

(2.3) H\? = V h ijia x«, 



\a\=k 



2 



(2.4) \H^)f = Y J Y J \ K 

ij \a\=k 

and for e > 0, set x = ey and define 

(2-5) H\f{y) = H^{ey). 

We will make extensive use of the following standard rescaling argument. Let {Ei}°2 =1 be a given se- 

pj- 1 

quence of positive numbers converging to zero. Define Mi by M- 2 = e, and in normal coordinates 

pj-i 

put y = M { 2 x = e i x and 

Vi {y) = M~ l Ui(x) = M^UiiM'^y) = ef~ x u^y) 



-i 



for y < crAI i 2 = e i 1 <r, where \x\ < a belongs to the domain of definition of the normal coordinates. 
Then v; satisfies 



(2.6) 



L h Vl + Kf~ 5i «f 1 = 0, for \y\ < aM^ 2 " . 

Bg,Vi = d v -Vi + ^K^Vi = 0, on dM, 



where U{y) = / 4 (M. * y) = f t (e iy ), fa) M (y) = (g^k^ \) = (giM^y) (see [251132]). 

We recall some standard definitions (see [251 [23]). Consider a sequence {ui} of solutions of (|2.ip . 
A point x € M is called a blow-up point for {ui} if Ui(xi) — > oo for some Xj — > x. 

Definition 2.1. A point i 6 Mis called an isolated blow-up point for {itj} if there exists a sequence 
{xi} C M, Xi — > x, where each X{ is a local maximum for Ui and 
1) uAxj) — > oo as i — >• oo, 
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2) Ui(x) < C dist 5i (x, Xi) Pi- 1 for x £ B a {xi) and some constants a, C > 0. 

Notice that the definition of isolated blow-up points is the same as for the boundaryless case (|23j). 

Remark 2.2. If we change the metric by a uniformly bounded conformal factor <j) > such that 
4>(xi) = 1 and V(j){xi) = 0, then isolated blow-up points are preserved. 

Definition 2.3. (|23j) Let {ui} and {xi} be as in definition 12, 1[ x^ — > x is an isolated simple blow-up 
point if for some p E (0,cr) and C > 1, where a comes from the definition of isolated blow up point, 
the function 

2 

2 r p i _1 /" 

■Ui(r) = rPi^uAr) = — - — — - — — — — — / u(z)dS(z) 
V1 y> vol 9i (MndB r ( Xi )) J MndBr(xi) yj y) 

Pi - 1 

satisfies, for large i, < for r such that CM i 2 < r < p. 

Observe that if x is an interior point then this definition agrees with the standard one (compare 
with [32]). 

2 — n 

Throughout the paper we let U : W 1 — >■ R be the function U(y) = (1 + \y\ 2 )^~ . U is known as 
the "standard bubble". From [25] we have the following. 



(2.7) A5 £ + n(n + 2)U^z e = c(n) £ didjH^U 



Definition 2.4. Let z £ be the solution of 

n-4 

/I/ . I j( fc ), 

fc=4 

constructed in [25j . It is implicitly assumed that z e = if n = 3,4,5. 
We recall estimate (4.4) of [25] 

n— 4 

(2.8) |^z £ (y)| < C J] ^e |a| |V Q |(l + |y|)l a l+2— 1/31 

H=4 £fc 

which implies 

(2.9) \d?z £ (y)\ < C{1 + \y\) 2 - n -W, for \y\ < ae~ l . 

The role of z £ is to provide a sharp correction term for the usual approximation of the (rescaled) 
solutions u by U around a blow-up point. z £ was introduced in the context of manifolds without 
boundary, and one of the main challenges in our paper is to establish that the same z E can be used 
in our setting. In other words, we need to show that z e satisfies a natural boundary condition. In 
order to accomplish this, we use one of the key results of the paper, theorem 13.41 to show that the 

~ (k) 

umbilicity of the boundary implies severe constraints on the behavior of the polynomials on the 

~ (k) 

boundary. Then we use the explicit construction of z £ in terms of to show that it satisfies the 
desired boundary condition. 

Notation and terminology used throughout the text: 

®d 

Pi- 1 

(ii) If X{ — > x is an isolated blow-up point, we denote M« = Ui(xi) and e~ = M i 2 . 

(iii) x' denotes the first n — 1 coordinate functions. 

(iv) We use iV to denote an integer that is arbitrarily large, coming typically from properties of 
conformal normal coordinates, such as det(g) = 1 + 0(r N ). 
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(v) Let Q be an open connected set that intersects dM. We then set d'Q, = 0,(1 dM and d + Q = 

dn\&n. 

(vi) In a coordinate system near the boundary, define -B^(O) = {x G B a (0) \ x n > G(x)} for 
some real valued function G, then denote d'B^(0) = {(x',G(x))}, d + B^(0) = dB^{0)\d'G^(0) (see 
section [3] and corollary I3.8P . 

(vii) We will always assume that the blow-up points x lie on the boundary dM, since theorems 
11.11 and II. 21 would otherwise follow from |25| (see section [T2|) . 

(viii) We will switch back and forth between problems (|2.ip and (|2.6p . referring to them as "x- 
coordinates" and "y-coordinates" . 

(ix) If xo € dM, then B a (xo) is a ball of radius a and center xo, i.e., 

B a (xo) = {x £ M \ dist(x,xo) < a}. 

Notice that B a (xo) will usually look more like a half-ball rather than like a full ball, but we will not 
denote it by B^(xq), reserving the latter for balls which explicitly satisfy the condition x n > 0. 

(x) For T > we define M" T = {x € R n \ x n > —T} and = {x € R n \ x n > 0}. 

3. Estimates near the boundary and boundary conformal normal coordinates 

In this section we will derive estimates for the second fundamental form, mean curvature, etc., in 
terms of the umbilicity tensor. Then we will use these estimates to modify the standard conformal 
normal coordinates construction in order to obtain conformal normal coordinates at the boundary 
with zero mean curvature. 

We first recall a result of Escobar. 

Proposition 3.1. (Conformal normal coordinates at the boundary \17\ ) Assume dM is umbilic and 
let xq € dM. For any N > there exists a metric g conformal to g such that, in normal coordinates 
for g centered at xq 

det(g) = l + 0(r N ), 

where r = \x\. If N > 5 then R s = 0{r 2 ), and AR- g (0) = -||Wg| 2 (0), and k s = 0(r 2 ). Here W s is 
the Weyl tensor. 

Now take conformal normal coordinates at xo € dM. We choose the coordinates so that di, 
i < i < n — 1, are tangent to dM and d n is normal (pointing inward) to dM at xo = 0. Near xq 
the boundary dM may be expressed as a graph x n = F(x'), where x' = (x\ . . . ,x n ~ 1 ) and since 
normal coordinates are defined up to a rotation we can assume that ^(0) = Vi ? (0) = and the 
tangent plane at is the "horizontal" hyperplane {x n = 0}. Then a basis for the tangent space 
T x dM is given by the vectors Xi = di + F^d n , 1 < % < n — 1. The normal is given as a covector by 
{vg)n = —1, {vg)i = 1 < i < n — 1, or as a vector by 

n-l 

(3-1) (v 9 y = 9 ij (u g )j = -g m ■ >>"/••,. 

i=i 

If g = e h then we may write 

{u g y = -8 in + h m + F >{ + 0(\h\ 2 + \h\\VF\). 
Define the second fundamental form by 

(3.2) By = K(Xi,Xj) = g{V Xi v g ,Xj), 



s 
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then 

(3.3) mj = giViVg, dj) + F ti g(V n Vg,dj) + Fjg(ViU g , d n ) + F ti Fjg(y n u g , d n ) 

= + F-ij - F ti T l nj (vg)i - FjV\ n (v g )i - FjF^T l nn (v g )i 

= \{-d n hij + dih nj + djhni) + Fij + 0(\h\\Vh\ + \VF\\Vh\). 
The mean curvature is given by 

n— 1 - 

(3.4) k = g~ 1 {X i ,X j )K i j = AF + £ 9 i h m + ^d n h nn 

i=i 

+ OQh\\Vh\ + |VF||V/i| + \h\\V 2 F\ + |VF| 2 |V 2 F| + \x\ N ), 
where we have used Y^=i ha = 0{\ x \ )■ Finally the umbilicity tensor is given by 

(3.5) T{Xi,Xj) = T l3 = niXuXj) - -^ K g{X i ,X j ). 

Notice that these quantities differ from the usual ones by a multiple of \v g \ (since v g is not necessarily 
a unit vector). As we show below (see proposition 13.51 and corollary I3.10p . this will immediately 
yield estimates for the standard (i.e., defined with respect to a unit vector) mean curvature, second 
fundamental form and umbilicity tensor, and it will suffice for our purposes. In fact, we will express 
all desired quantities in terms of Tij, and the umbilicity of the boundary implies that Tij defined 
with respect to (|3.ip vanishes as well. We remark also that our definition of the mean curvature in 
this section differs from the standard one by a multiple of (n — 1) . However, in all other sections 
of the paper we adopt the standard convention, unless otherwise specified. 

The next theorem will be our main tool to produce estimates. Although its proof is long, the 
idea behind it is quite simple: from properties of conformal normal coordinates we can derive several 
identities involving geometric quantities and the functions hij. We restrict the obtained expressions 
to their Taylor polynomials, and successively solve these equations for one quantity in terms of the 
others, until we express all quantities in terms of the umbilicity tensor and an error. 

Remark 3.2. It should be noted that in (j3.3[) and (|3.4|) . as well as in the proof below, the expression 
|/i||V/i| appearing in the error only includes terms of the form |/i||9j/i n j|, |/i||9 n /ijj|, |<9 n /i nn | or 
\h ni \\Vh\. 

Remark 3.3. Since we will eventually restrict all expressions to their Taylor polynomials in theorem 
13.41 and is arbitrarily large, we will ignore the 0(\x\ ) contributions. 

Theorem 3.4. Take conformal normal coordinates at xq G dM as described above and choose a large 
integer N. Then there exists a constant C , depending only on N such that for any e > sufficiently 
small: 

N N n-1 

\a\=2 \ a \=2i,j=l 
N N-2 n-1 

£ |AF a | £ H < C £ \T^a\e^ +2 

\a\=2 [ a |=Qij'=l 
N n-1 N n-1 

H=0i,j=l |a|=0*J=l 
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N N-2 n-1 

\a\=2 |a|=0«J=l 
N n-1 AT-1 n-1 

EEi^i^'^E y,\ t ^ h+1 

|a|=2i=l \ a \=H,j=l 
N n-1 N n-1 

E E Idnhij,*^ 1 < C E l%«l e ' a| 
a|=li,j=l |a[=lij=l 

AT N n-1 

E l<9n/wk H < c E E l^>k |a| 

|«|=1 | a |=lij=l 

where a denotes partial derivatives in the variables x x n_1 evaluated at the origin, and F is the 
local representation of the boundary as a graph as explained at the beginning of this section. Moreover 
k(0) = |V/c|(0) = F(0) = |VF|(0) = AF(0) = and |VAFj(0) < C£y |V^|(0). 

Proof. We first record several useful calculations. When repeated indices i or j appear this signifies 
summation from 1 to n — 1. Using familiar properties of conformal normal coordinates and (|3.3[) we 

have 

(3.6) = i[-a„(x'^) + ,% Kj - S> lKi + + xV« + OdxIIMIVfcl + M|VF||V*|) 
= ^x l d;/i nj +x i d i F d + 0{\x n \\Vh\ + |x||/i||V/i| + |x||VF[|V/i|). 

Furthermore 

(3.7) x x^ Kij = x x^ Fij —x S^h n j -\- — x di(x^h n j) 

+ 0(|x[|x n ||V/i| + |x| 2 |/i||V/i| + |x| 2 |VF||V/i|) 

= iV'f^ + ^x n /i nn - -x^dihrm + 0(\x\\x n \\Vh\ + |x| 2 |/i||V/i| + [x| 2 |VF||V/i|) 
= xV>y + 0(|x n ||/i| + |x||x n ||V/i| + |x| 2 |/i||V/i| + |x| 2 |VF||V/i|), 

and 

(3.8) x l x j Kg(Xi,Xj) = x l x j K(g i: j + F^g jn + + F^Fjg nn ) 

= \x'\ 2 K + 0{\x\ 2 \h\\\ K \ + |x| 2 ||VF| 2 |k|). 

Recalling the definition of the umbilicity tensor together with (|3,7p and (|3.8p yields 

(3.9) xVFy = -iy|i'| 2 K + iV7ij +0(|x"||/i| + |x||x n |[V/i| + |x| 2 |/i||V/i| 

+ |x| 2 |VF||V/i| + |x| 2 |/i||k| + |x| 2 |VF| 2 |k|). 

Moreover since 

x i ng{X i ,X j ) = x j K + 0(\x\\h\\K\ + |x||VF| 2 |k|), 
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we find that (using (13. 6p ) 

(3.10) -x^iKj+x^iFj 

= —-jX^K + x^j + 0(\x n \\Vh\ + \x\\h\\Vh\ + |x||VF||V/i| + \x\\h\\K,\ + |x||VF| 2 |k|). 
Eliminating k from (j3.9|) and (|3.10p produces 

(3.11) x^iFj + -x^ihnj - x% = \x'\- 2 x j x i x\F u - T u ) + Slj, 
where throughout this proof fL denotes error which satisfies 

(3.12) = 0(\x 1l \\Vh\ + |x||/i||V/t| + |x||VF||V/i[ + \x\\h\\K\ + \x\\VF\ 2 \k\ + |z , |- 1 |z n ||^|) 

= 0(\x n \\Vh\ + |x||/i||V/t| + |x||VF||V/i| + |x||/i||V 2 F| + |x||VF| 2 |V 2 F| + |x / | -1 |s n ||/i|). 
Upon restricting attention to Taylor polynomials (|3.1ip simplifies to 

(3.13) {k - l)\x'\ 2 Ff - k{k - l)x j F^ 

= -\{k - l)\x'\ 2 h^ + [xfx'T^ - xix^T^ + \ x '\ 2 nf- l) 

where h^j ^ denotes the (k — l)-degree Taylor polynomial in the variables x 1 , . . . , x n_1 and similarly 

for F^ k \ T^j 2 . We note that h^- ^ is not the full Taylor polynomial in all the variables x 1 , . . . , x n 
but rather just the portion involving the first n—1 coordinates, and the remainder involving x n is 
relegated to the error term. Now apply dj to f)3. 13|) and sum over j to find an equation for F^ k \ 

(3.14) \x'\ 2 AF^ +k(3-n- k)F^ = -^\ x fdjh^ 

+(k - l^djdxfx'T^ - xVx'lf^) + d j (\a?\ 2 nf- 1) ) 

where we used x l hij = —x n h n j to absorb this term in the error. Differentiate (|3.10p with respect to 
x 3 and sum over j to find 

^djh n j + ^x l didjh n j + AF + x l diAF = x % d.fT^ + k H —x J djK + djVtj 

where we used that the error term in (|3.10p has the form Qj. Then 

(3.15) \{k - 1)8^ + (k — 1)AFW = Jdfl§-*> + n + k _-\ ^2) + d . n (k-D_ 
On the other hand (|3.9p gives 

(3.16) k(k - 1)F« - xV'ri*" 25 = — |x'| V fc ~ 2 ) + x^V. 

Tl 1 

Therefore using (|3.14p and (|3.16p in (|3.15p produces 

(3.17) 8>r> = -eV* 1 ^ + w^ 1 ' + 

Let £>™ _1 denote the unit ball with respect to x l ,...,x n ~ l , and let 4> G C 00 ^ -2 ). Extend 4> 
radially so that it is defined on B™~ 1 \{0} and d r <p = on S™~ 2 , where r = \x'\. Notice that even 
though (j> is not defined at the origin, we can still integrate by parts against functions which vanish 
at zero, and so in particular against homogeneous polynomials. 
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Let (j) be as above. From (|3.10j) we have 

(3.18) (k - 1)8^) = — !— j-a^K^ -2 ) - + x^ 2) + nt 1] - 

Multiply (|3.18p by di<fi, sum over i and integrate by parts to get 

-(fc-1) / (f>AF^ + {k-l) f (j>u%F^ = — / Mn + k - 3)^-^ 

Jb^' 1 Js™~ 2 n-l Jb™- 1 

+^-r [ ^xV*- 2 ) + k —± I Mt^-V _ til [ 



77, — 1 J s n-2 z J B n ~ * J S n ~ 



Bp 1 JS™" 2 JBp 1 JS™" 2 



Integrating in polar coordinates produces 

As™" 2 Js™ 



(Jfe - 1) / r n+k ~ A [ <f,AF^ + (k - 1) / ^diF^ 
Jo Js™- 2 Js™- 2 

1 rfc+ „.„4 /" 0(n + fc _ 3)K (fc- 2 ) + 1 /" ^ 
™ - 1 Jo Js?- 2 n ~ l Js\ l ~ 2 



V i X i K<- k - 2 '* 



+ k-i r 1 rn+k _ 4 r x) _k-i f ^ v _ r 1 r „ +fc _ 4 . ^ 



y^r 2 jo isr 2 ^r 2 

which implies (notice that the mean curvature terms cancel out) 

(3.19, afw - ^^±^f<" = -i^- 11 + i 1 ^^-" + ^T^8^.r 2) 

n + k-3iV ( fc _2) , o„(fc-i) , ^ (*-i) 

where we have used that 4> is an arbitrary smooth function on S™ -2 and homogeneous polynomials 
are determined by their values on the sphere. Using (|3.14p and (|3,17p in (|3. 19j) we find that 

2 

(jfe-l)(n + ife-3) 
Similarly, multiplying f)3. 13j) by dj<f> and integrating by parts yields 
(3.21) 

- ^^xW^+x^f^ + Ixfd^. 
Solving for AF^ + hdjh^f 1 ^ in (jHJH) and using it along with QHJ2CTD in ([3T2T]) we obtain 



(3.20) X ^ = - (h _, v ' h _^ T^Tt 2) + + l^-n^. 



X 



Notice that the denominator of the first term on the right hand side is never zero since k > 2. 
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From (|3.9h we have 

(3.23) k\x'\- 2 F^ = 7 1- + J—\ X '\-^T^ + 

(n-l)(k-l) k — 1 1 3 

Using (^22]) in ([333]) yields 

where c(n, fc) is a numerical factor depending on n and k only. Let ft be the set of homogeneous 
polynomials that can be estimated in terms of the umbilicity tensor and an error (of the same degree) 
in Qj. Then (1337]) . (13301) . (13321) and (1331) give that d^'^ , x*/^" 1 ^ F^ k \ k^'^ G ft. From 
(13"T9]) it then follows that AF^ G ft as well. From (13TT3j) and S K we get G ft, and from 

dS3D and AF&\ dih%~ l) G 72. it follows that (d n /i n „)^- 2 ) G ft. Using ((33]) along with ^ k ~^ G ft 

we get Kfr 2) e ft and from this, (33]), /iJJ - ^ G ft and M fc ) G ft we find that (fl^/iy) ( *~ 2) € ^- 
The inequalities of theorem 13.41 now follow with the help of remark 13.21 

By our construction of F and properties of conformal normal coordinates we have k(0) = | Vk| (0) = 
F{0) = \X7F\(0). Hence in order to finish the theorem we only have to show that AF(0) = and 
|VAF|(0)<C£ v -V|T^|(0). 

Using the definition of K{j, and recalling that v n = — 1 and Uj = Fj, 1 < j < n — 1, we obtain 

(3.25) K ij = F ,ij ~ T ij F ,k + T™- - F tj T k n F jk + -FjT™„ 

— F T k Fi 4- FT n — FF T k Fu+FFT n 

where we have used F n = and X)fc=i Ffj-^fc = ^'■^A smce F does not depend on x n . Evaluating 
([335]) at and using r*-(0) = = |VF|(0), we have Ky(0) = ify(0). Taking a trace produces 
AF(0) = k(0) = 0. Finally notice that (f3TT2l) gives 

fii = 0(|x ra ||x| + |x| 4 + |x| 3 + |x'r 1 |x n |k| 2 )> 
so we can compute directly from f)3.22[) to find |VAF(0)| < C^ij |V2#(0)|, finishing the proof. □ 

Now with the help of theorem 13.41 we improve the properties of conformal normal coordinates 
at the boundary by showing that we can also require zero mean curvature. We call these coordi- 
nates boundary conformal normal coordinates to avoid confusion with the usual conformal normal 
coordinates at a point on the boundary. 

Proposition 3.5. (Boundary conformal normal coordinates) Let (M, go) be a Riemannian manifold 
with umbilic boundary and xq G dM . Fix an integer N > 5. Then there exists a metric g conformal 
to go such that, in g-normal coordinates centered at xq: (i) detg = 1 + 0(r N ), (ii) R g = 0(r 2 ), (Hi) 
A g R g (0) = — ^\W g \ 2 (0) and (iv) Kg = near xq, where r = \x\. 

Proof. Using conformal normal coordinates at xq we obtain a metric g which satisfies properties 
(i)-(iii) in a ball B a (0). Our task is to show that we can perform a further conformal change in the 
metric in order to obtain property (iv) while maintaining (i)-(iii). 

We write all quantities as explained above (see equation (j3.ip and what follows); in particular we 
denote by n g the mean curvature defined as in (j3,4p . and by n g the mean curvature defined in the 
usual way, i.e., with respect to a unit vector. 

If g = e 2 ^ g then K g = e~f(K g + Jjj-). We will choose / appropriately. 
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Because the boundary is umbilic, Ty vanishes identically and therefore theorem 13.41 gives K g>a (0) = 

for \a\ = 0, . . . , N, where a denotes derivatives with respect to x , . . . , x n ~ l . In other words, we 
obtain that K g = 0(\x'\ N ), from which it follows that n g = 0(\x'\ N ) as well. Now we choose an 
extension of K g to H g in a neighborhood of xo, with H g satisfying H g = 0(1x1^) and = . Such 
an extension is possible because K g = 0(\x'\ N ). 

Now pick a smooth function / such that J^f- = — 1 near xq and put / = fn g . With this choice of 
/ we then have K g = in a neighborhood of xq. 

By construction we have / = 0(1x1^), and so we obtain the desired result as the remaining 
properties all follow from det(g) = 1 + 0(r N ) (after choosing a smooth extension of / to the whole 
ofM). □ 

Remark 3.6. We stress a point made in the introduction. The so-called conformal Fermi coordinates 
|33j have been used with great success in the study of the Yamabe problem for manifolds with 
boundary (see references mentioned in the introduction). This expresses the fact that cylindrical 
coordinates generally work better than spherical ones for Neumann-type of problems. However, a 
critical part of the compactness result of Khuri, Marques and Schoen [25] for boundaryless manifolds 
is the proof of the positivity of a quadratic form on Taylor polynomials of the scalar curvature 
which naturally arises in the problem. Their proof makes substantial use of the the radial symmetry 
coming from normal coordinates and we would like to preserve as much as possible of that original 
argument. Boundary conformal normal coordinates preserve the radial symmetry while displaying 
features similar to the good properties of Fermi coordinates, as it is shown below. 

Boundary conformal normal coordinates have the following useful property. 

Corollary 3.7. In boundary conformal normal coordinates centered at xq € dM the boundary is 
given by x n = 0. Moreover, gi n (x',0) = 0(\x'\ N ), 1 < i < n — 1. 

Proof. Since the boundary is umbilic and n g = 0, it is also totally geodesic (i.e. Kij = 0) for the 
metric g, and therefore the boundary is given by x n = in normal coordinates. The second statement 
follows from theorem 13.41 as g = e . □ 

Now we want to extend the previous results for the case of interior points. Assume that xq € M 
is an interior point which is sufficiently close to dM, and take conformal normal coordinates at xq. 
Denote by Jo € dM the closet point to xo- We can still write the boundary as a graph x n = F(x'), 
and since normal coordinates are defined up to a rotation we can assume that -F(O) = — |xo| where 
|x | = dist(xo,xo) ( so that xo = (0, . . . , 0, — |xo|)) and the tangent plane is horizontal there, so 
|V-F(0)[ = 0. Moreover, by the Gauss lemma we also have = g nn <9 n | £ . 

If we "translate the boundary", i.e., define 

G(x') = F(x') + |x | 

we have G(0) = |VG(0)| = and d a G = d a F. Set B% = {x € B a (0)\x n > G{x')}. Notice then that 
a basis for the tangent space at a point on d'B^ = {(x',G(x'))} is X{ = d{ + Gjd n = di + F^d n , 

1 < i < n — 1 (since we are simply translating the boundary). We can then consider all geometric 
quantities induced on the boundary d'B^. In this situation, theorem 13.41 holds with G replacing F 
and all quantities being defined with respect to the boundary d'B^, except for the conclusions that 
depend on k = 0(r 2 ), since the boundary d'B^ need not be umbilic. We state this as a corollary. 
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Corollary 3.8. Let xq € M , take conformal normal coordinates at xq and assume that xq is suf- 
ficiently close to dM as to have dM n B a (0) ^ 0, where B a (0) is the domain of definition of the 
conformal normal coordinates. Let F = F(x') be the local representation of the boundary as explained 
above. Define G(x') = F(x') + \x \, B% = {x G B a (0)\x n > G{x')}, and d'B% = {(x 1 ,G(x'))}. Then 
there exists a constant C , depending only on N , such that for any e > sufficiently small: 



N N n-1 

E m^ 1 <cE E 

H=2 |q|=2«J=1 
N N-2 n-1 

y: iaguism < c y: e 

\a\=2 | a |=0i,J=l 
JV n-1 TV n-1 

Y E ^ E E 

|a|=0i,j=l |a|=0ij=l 
N N-2 n-1 

E \ G ^ H ^ E E 

|a|=2 | a |=0i,i=l 
N n-1 N-1 n-1 

|a|=2j=l \a\=li,j=l 
N n-1 N n-1 

E E i^'> |Q| < C E E 

|a|=li,i=l \a\=H,j=l 
N N n-1 

E \dnhnn,a\e H < C E E 
|a|=l |a|=l ij=l 



- 1 | c 



,M+2 



|a| 



r |a|+2 



r |o|+l 



rl«l 



rl«l 



where k, Kij and Ty are respectively the mean curvature, second fundamental form and umbilicity 
tensor of &B^ , all defined with respect to the outer normal 

n-1 
i=l 

(which is not necessarily a unit normal) and a denotes partial derivatives in the variables x , . . . , x™" 1 
evaluated at the origin. Moreover G(0) = |VG|(0) = and |VAG|(0) < CJ2ij |V2^|(0). 

As before, estimates on quantities defined with respect to u g , with v g not necessarily a unit vector, 
will suffice for our purposes. 

Because d a G = d a F, estimates for G from corollary 13,81 translate into estimates for F. 

Corollary 3.9. Let xq £ M and xq € dM be the closest point to xq- Take conformal normal 
coordinates at xq, choose a large integer N and let F be the local representation of the boundary as 
a graph. Then there exists a constant C , depending only on N such that for any e > sufficiently 
small: 

N N-2 n-1 

£ |F a | £ H < C E E l^ H+2 

H=2 | a |=0M=l 
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£ |AF a | £ H < c £ £ \n 



e 



a\+2 




where is the umbilicity tensor of the boundary d'B^ as in corollary \3.8l and a denotes partial 
derivatives in the variables x , . . . , x n_1 evaluated at the origin. Moreover IVFKO) = 0, |VAF|(0) < 



The following corollary will finish the treatment of interior points in this section. 

Corollary 3.10. (Boundary conformal normal coordinates for an interior point) Let (M,go) be a 
Riemannian manifold with umbilic boundary and xq G M . Fix an integer N > 5. If xq is sufficiently 
close to dM , then there exists a metric g conformal to go such that, in g-normal coordinates centered 
atx : (i) det<? = l + 0(r N ), (ii) R- g = 0(r 2 ), (Hi) AgR g (0) = -||Wg| 2 (0) and (iv) K g = near x , 
where r = \x\ and xq € dM is such that dist 9o (a;o> ^o) = dist 9o (xo, dM). 

Proof. Let xq G dM be the closet point to xq. Denote by {5*} conformal normal coordinates 
centered at xo, {x 1 } conformal normal coordinates centered at xo, k the mean curvature of dM in 
{x*}-coordinates, k the mean curvature of dM in {x*}-coordinates. When xq — > xq we have x l —} x l , 
and d a K(xo) — > 9qk(xo) where a denotes partial derivatives with respect to x l ,...,x n ~ l and a 
denotes partial derivatives with respect to to x , . . . , ac n—1 . 

By theorem 13.41 we have that dan(xo) = for |d| < N since the boundary is umbilic. Because 
d a n(xo) — > dan(x~o) as xq — > xq, if xo is sufficiently close to xq we can choose an extension of k to 
B a {xo) (with a small) which is 0(\x — xq\ n ). The rest of the argument now is similar to the proof 
of proposition 13.51 □ 

We finish this section with several remarks. 

Remark 3.11. One of the key ingredients of our proof is to show that the blow-up sequence Xi 
lies on the boundary (possibly after passing to a subsequence, see section [7]). Before showing that, 
however, we have to deal with both the case of a blow-up sequence belonging to the boundary and the 
case of a blow-up sequence belonging to the interior of the manifold. It will therefore be implicitly 
understood that when Xi £ M, all quantities At, kij and are for the boundary d'B^, as described 
above, i.e., we will drop ~ from the interior quantities for the sake of notation. F, however, will 
always be the representation of dM as a graph unless stated otherwise. 

Remark 3.12. Suppose that xq € dM or that it is sufficiently close to the boundary, and in 
boundary conformal normal coordinates centered at xq consider x = (0,x n ). If we translate the 
boundary by \x n \ instead of \xq\, 



we can, for each \x n \, consider geometric quantities induced on the boundary d'B^ as before. In 
another words, we have a foliation of a small neighborhood of the boundary by copies of dM. 
In particular, we can then think of as defined in a neighborhood of dM, allowing us to take 
derivatives with respect to x n , Taylor expand Tij in the x n direction, etc. 

Remark 3.13. Since boundary conformal normal coordinates are a special case of conformal normal 
coordinates, the results of this section stated for conformal normal coordinates, in particular theorem 
13.41 are still valid if we choose boundary conformal normal coordinates instead. 



CE«|V3ij(0)| and F(0) = -\x 



G(x') = F(x') + \x 
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4. Higher order estimates 

Our next goal is to extend the results of theorem 13.41 to higher order derivatives of h in the 
normal direction. Throughout this section we will work with boundary conformal normal coordinates 
centered at a point on the boundary; all definitions are as in section [3l We will use Greek letters 
to denote indices running up to n, Latin letters to denote indices running up to n — 1, and x' to 
denote the first n — 1 coordinates. Notice that in light of corollary 13.71 we have that the boundary 
is given by x n = 0, and as in section [3l by a "translation" we can consider quantities defined on the 
neighborhood of the boundary, so that 

(4.1) /- = -g nT d T , 

(4.2) gni\ dM =9m(x' 1 0)=O(\x'\ N ), 





(4-3) 

" V 9 



dv n 



g nn (x',0)d n + Y,O(W\ N )de- 



dM e=1 
Theorem 4.1. In boundary conformal normal coordinates at a point on the boundary, 

l nn\QM 



(4.4) h nn \ dM = h nn (x',0) = O(\x'\ N ) 



Proof. We will compute 'S/iV n in two different ways. First, 

(4.5) Viu n = -d t g nn - g™T? n - g nl T^. 

Notice that (|3.3[) becomes in our coordinates T^(x',0) = Kij(x',0) = 0, and hence (|4.5p gives 

(4.6) ViV n \ dM = -d ig nn - \{g nn ?d ignn - \g nn g n \d ignl + d n9li - d l9m ) 

= -d l9 nn - \(g nn ) 2 d l9 nn + OdxT- 1 ), 

where we used theorem 13.41 In order to simplify notation, here and in the rest of the section we use 
the following convention. When an equality is restricted to the boundary we write -\ dM or (-)(x',0) 
on one side of the equation, and it is implicitly understood that the remaining quantities on the 
other side are restricted as well. 

Now differentiate g nn g nn + g nl ' g n i = 1 with respect to i to obtain g nn dig nn = - g nn dig nn - di(g nl g ni ) 
and so (|4.6p becomes 

(4.7) V^U = -d l9 nn + l -g rm gnnd ig nn - d^gm) + 0{\x'\ 2N - 1 ) 

= -d i9 nn + \(l - g nl g n i)d t g nn - W'snz) + Odx'l^" 1 ) 

= -ld t g nn + 0(\x'\ 2N - 1 ), 

where we used theorem 13.41 again. Combining (|4.6p and (|4.7p gives 

d ig nn + {g nn ) 2 d ignn = 0{\x'\ 2N - 1 ). 

Using g = e h this becomes 

(4.8) -(1 - (g nn ) 2 )dihnn + diO nn (-h) + (g nn f diO nn (h) = Odxf^" 1 ) 
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where 



But 



and 



o^)=f: ihe) 



O n n(-h) = -(h nn h nn + h n ih n i) + ' (21)"" ~~ 5^ ^(21 + 1)" 



I — (h 2e ) nn °° (h? u 

O nn (h) = -{h nn h nn + hnihni) + ^2 (npT + E 



2^ 

nn 



2 v-....--.«. ■ ■ ^ (2^; (2^ + 1)!' 

and since by theorem El /i n j(a/,0)/i Bl (sc / ,0) = 0(\x'\ 2N ), gSJ) becomes 

(4.9) (-1 + G?"") 2 + h nn + (g nn ) 2 Kn)d l h nn + (i + GT) 2 )^ V 

+(i - (ff -) 2 )^ E = od^'i 2 ^ 1 )- 

We will now show inductively that h nn = 0{\x'\ k ) implies h nn = 0(\x'\ 3k ), k < N. Notice that we 
already know that h nn = 0(|x'| 2 ). Also, as before, the terms h n i appearing in (h 2i ) nn , i > 2, and 
(h 2e+1 ) nn , £ > 1, can be estimated by theorem 13.41 and hence they can be absorbed in the error; in 
other words we can replace (h £ ) nn by (h nn ) e up to and error 0(|x'| 27V_1 ) (notice that due to the rules 
of multiplication of matrices, the terms h n i appearing in (h 2i ) nn , £ > 2, (h 2i+1 ) nn , £ > 1, or in the 
expansion of g nn must be multiplied by another h n i and hence such errors are of the same order of 
the right hand side of (|4.9p ). 

Since g nn > C > near the origin, if h nn = 0(|x'| fe ) then 

-1 + {g nn f = -(1 + g nn ){\ - g nn ) = -(1 + g nn )0{\x'\ k ) = 0{\x'\ k ) 



and also 



Kn + (g nn ) 2 Kn = 0(\x« k 



SO 

(4.10) -1 + (g nn f + Kn + {g nn ?Kn = 0(\xf). 

But 

(4.H) (i + { 9 n 2 )d* £ = (i + (g nn m f; %2 + ow™- 1 ) 

1=2 1=2 

= 0(\x'\ 4k - 1 ) + 0(\x'\ 2N - 1 ) 

and 

(4.i2) (i - {g-m f; = (i + «rxi - <™ f; + 0(\xT-l) 

1=2 ^ >' 1=2 ^ >' 

= (1 + <7 nn )0(|x'| fe )0(|x / | 3fc - 1 ) + Odx'l 2 ^ 1 ) 
= 0(|x'| 4fc - 1 ) + 0(|x'| 2W - 1 ). 
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Therefore P~gj) - P~T2j) give 



d-h - - fl I (a nn ) 2 )d V ^ )nn 

i nn ~ (-1 + GrO 2 + c + (r) 2 M ( + 15 j M h W 



1=2 
00 



+ I fi _ (n nn ) 2 )d- V ( /i2€+1 ) 



(-1 + (S n ") 2 + ^nn + (5 nn ) 2 /inn) ' ^ ' ' ^ W + W 

= Odx'l 3 ^ 1 ) + 0{\x'\ 2N - 1 - k ) : 

provided that (-1 + (g nn ) 2 + h nn + (g nn ) 2 h nn )(x' , 0) is not zero and k < N. Since /i(0) = we 
conclude that h nn (x') = 0(\x'\ 3k ). Repeating the argument we obtain the result. 

Now we have to show that the result is still true if (—1 + (g nn ) 2 + h nn + (g nn ) 2 h nn )(x' , 0) vanishes 
or is 0(\x'\ N ), and it is enough to consider this latter case. So suppose that (—1 + (g nn ) 2 + h nn + 
(g nn ) 2 Kn)(x',0) = 0(\x'\ N ). Multiply it by g nn and use 1 = g nn g nn + g nl g nl = g nn g nn + 0(\x'\ 2N ) 
to get 

-9nn + + (9nn + Q^Kn = 0(\x'\ N ), 

which implies -O nn (h) + O nn (h) + (O nn (h) + O nn (-h))h nn = 0(\x'\ N ) and therefore 

2 V ( h2t+1 )nn 2h {h 2i ) nn _ Q,y,N\ 

{2i+l)] +^nn^ m ~0(\X\ ). 

As before we can ignore contributions from h n i and replace (h l ) nn by (h nn Y, which gives 

- oo /, 2^+ 1 00 1 ^ \2£ 

-\(h nn f - 2 £ ^f- - + ((/.nn) 2 + 2 £ ^f) ^ = 0(|,r )• 

This gives (h nn ) 3 = 0((h nn f) + 0(\x'\ N ). Since ft = 0(\x\ 2 ) we obtain (h nn f = Od^'l 10 ). But then 

(hnnf = 0((h nn f) + 0(\x'\ N ) = 0((h nn f(h nn f) + 0(\x'\ N ) 

= 0(\x'\ l \h nn ) 2 ) + 0(\x'\ N ) = 0(\x'\ 1A ). 

Repeating the argument produces (h nn ) 3 = 0(\x'\ N ), which gives the result since N is as large as 
we want. □ 

Theorem 4.2. In boundary conformal normal coordinates centered at a point on the boundary we 
have 

(4.13) d 2 n h nj , al (Q) = 0, \a'\<N 

(4.14) d*h nn ,a'(0)=0, \a\ < N 

where a' denotes derivatives with respect to x 1 , . . . ,x ra_1 . In other words d 2 ,h n j\ [QM = 0(\x'\ N ) and 

dn,hnn\g M = 0(\x \ N ). 

Proof. Denote by h*f? the m th Taylor polynomial of h n i. Let <fi £ C°°(S'™ _1 ) and extend it radially 
similarly to what was done in theorem 13.41 (notice however that here we have the full, i.e., including 
x n , Taylor polynomial). Integration by parts yields 

2 d nh ^ = -2[ 4^0$ + [ y^-f yh^ 

B+ JB+ JS^- 1 Jb™- 1 
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where B + is the half unit ball and B n ~ l the unit ball in x' coordinates. Since h n i(x',0) = 0(\x'\ N ) 
by theorem 13.41 we obtain that the integral over _B n_1 vanishes. Integrating in polar coordinates as 
in theorem 13.41 shows that 

tfajtf* = - 2 m + n ~ l [ + (m + n - 1) [ <?x n h^\ 

-i m m + n Js™- 1 

And since this is true for any e C 00 ^- 1 ), we conclude 

4*nt& = -2^±^d n ^ + *(m + n - l)*h™ on S^. 

Using theorem 13.41 again, or, alternatively, choosing a non-zero test function such that d n (j) = on 
gn-2 _ Qgn-l^ j t follows that d n h^\x',0) = 0, from which we conclude 

(4.15) d n h nl (x',0)=O(\x'\ N ). 

Now with (|4.15p in hand, we repeat the integration by parts argument with d\h n i in place of d n h n \ 
and conclude (|4.13p . 

To obtain (|4.14p . argue similarly to the above, integrate d\h nn by parts and use theorem 13.41 to 
conclude d\h nn (x\ 0) = 0(\x'\ N ); then repeat the argument, expressing d^h nn in terms of d^h nn . □ 

Remark 4.3. Since theorem 13.41 gives d n h nn (x',0) = 0(\x'\ N ), using an argument similar to that of 
theorem 14.21 we can relate d n h nn (x', 0) and h nn (x',0), obtaining in this way an alternative proof of 
theorem 14.11 

5. Boundary condition for the correction term 

In this section we use the results of sections [3] and [J] to show that the correction term z £ (see 
definition I2.4p satisfies the correct boundary condition. The idea is to use the results of sections [3] 
and [4] to show that certain homogeneous polynomials that appear in the (explicit) construction of z e 
satisfy the boundary condition and so will z e itself. Throughout this section we work with boundary 
conformal normal coordinates centered at a point on the boundary. This section relies heavily on 
the appendix of [25J and we will often refer to it. 

Lemma (A. 6) of [25J gives the decomposition 



[^1 

( fc ) _ w« _ 



(5.i) 

9=1 

where satisfies <% W^- = 0. In [25] it is also computed that 

(5.2) = " 2q)(k -2q- l)(n + k — 2q - l)(n + k-2q- 2)|zf«- V-2 9 

H n — 1 

= Cn,k,q\ x \ 2q 2 Pk-2qi 

where Pk-2q is a harmonic polynomial of degree k — 2q. 
Since H^\x) = Yl\ a \=k h ij,aX a , we obtain 

^n^ijH— (x) — ^ ^ hij jCt O n ijX . 



Now we consider d n dijH\- 



\<x\=k 

(k) 



and identify the terms that do not necessarily vanish (recall that the 
boundary is given by x n = 0). 
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Consider the case i,j < n. In this case if 

^ ^ hij,a 9ij n X \ x n_Q ^ 0, 
\a\=k 

then the non-zero terms have a n = 1, i.e., we can write the multi-index a (of the non-zero terms) as 
a = (a 1 , 1). Hence the coefficients of the non- vanishing terms are all of the form 

(5.3) h ij)Ct = ^-d a h ij (0) = ^-d n h ij a ,(0), \a'\ = k - 1. 

Similarly if i = n and j < n then the coefficients of the non- vanishing terms are all of the form 

(5.4) h ijj0 , = ^dlh nj ^(0), \a'\ = k - 2, 

and if % = j = n the coefficients of the non- vanishing terms are all of the form 

(5.5) h ij a = —8%h nn a/ (0), \a'\ = k - 3, 

where a' is a multi-index with a' n = 0. Since k < n — 4, we have by theorems 13.41 and 14.21 that (|5.3p . 
(|5.4p and (|5,5p all vanish, and therefore 



(5.6) d n d ij H^\x',0) = 0. 



Combining (|5.ip and d^W^ = with (j5.2j) and (|5.6p gives 



C nAq \x'\ 2q - 2 d nPk „ 2q {x\Q) = 0, 

9=1 

and it then follows from usual decomposition theorems for homogeneous polynomials (see e.g. [6]) 
that each d n pk-2q(x' , 0) vanishes separately. 

To compute d n z e \^n-i it is enough to compute the derivative of Z((Hq k ^)ij) — the solution to 
(|2.7p with (Hg k ^)ij instead of J2k^-ij^- Such a solution takes the form (see [25] ) 

9+1 

Z({HW)ij) = c(n)a k _ 2q (l + |xj 2 )~t ]T T(k, q, j)\x\ 2j Pk-2q 

3=1 

where a k -2q and T(k,q,j) are numerical coefficients. Computing we find 

9+1 

d n Z((i?«)y)(x',0) = c(n)a fc _ 29 (l + |xf )"f ^r(fc, g, j)|xf ^ nPfe _ 2? (x', 0). 

i=i 

But we showed above that d n pk-2q(%', 0) = and hence d n Z((Hq )ij)(x', 0). Therefore, we have 
proven 

Proposition 5.1. Ta&e boundary conformal normal coordinates at a point on the boundary and let 
z e be as in definition \2.J\ Then it satisfies 

(5.7) 



' Az £ + n{n + 2)U—~z £ = c(n) YTklt didjH^U, in 



d n z £ = 0, on 



n—l 
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6. Basic convergence results 

Here we prove some basic convergence results. Most of the results are either known or modifications 
of similar results for manifolds without boundary. 

Lemma 6.1. Suppose x% — > x is an isolated blow-up point. Take normal coordinates at x%, and 
rescale coordinates to y- coordinates. Then in the limit i — > oo the boundary becomes a hyperplane. 

Proof. The metric gi is obtained from g by (i) a rescaling g\ i = Mf 1 ~ 1 gi = £^ 2 g and then (ii) by 

4 

the change of coordinates y = e x. If we write g\ t = M^'~ g^ in the standard form g\. = 4>™~ 2 gi 
_JL_ Pi- 1 

we have 4>™ 2 = M i 2 = e i , so transformation law ()A.4j) gives (rei^jy = £j ( K i)kj- The second 
fundamental form transforms as 

dx^ dx^ 

= £ l(. K u)kj(x) = £iE~ l (hi)kj(x) = ei(Ki) kj (x) 

when we change coordinates from x to y via x = e^y. Now the sequence (Ki)kj(x) is bounded because 
in x-coordinates the metrics converge in C k (k large), and therefore the second fundamental form 
goes to zero in y-coordinates. But gij{y) — > 5ij since 

9ij(y) = 9ij{£iy) = Sij +eiO{\y\ 2 ) 
and therefore in the limit the boundary is a hyperplane (see also |23j). □ 

Lemma 6.2. Let Xi — > x £ dM be an isolated blow-up point. There exists a constant C > such 

Pi- 1 

that for alii and all \y\ < o~M i 2 we have \vi(y)\ < C (where a comes from the definition of isolated 
blow-up point). 



Proof. The proof is similar to the first claim in proposition 4.3 of |32j and it uses the maximum 
principle, the Harnack inequality and the definition of isolated blow-up point. In fact, from the 
definition of vi and isolated blow-up points we have that 



3.1) 



'vi(0) = 1, Vwi(O) = 

.0 < Vi(y) < C\y\~^ for \y\ < aM~ = U. 



i', 



From these properties it follows that i>i(y) < C for 1 < \y\ < l{. Since Lg t Vi = —Kvf 1 < 0, using the 
maximum principle (corollary IA.4|) we have that there exists a constant C > such that for every i, 

min Vi(y) > min Vi{y) 

\y\<r \y\=r 

with < r < 1. Using the Harnack inequality (lemma |A.6|) we get 



so that 



m&xvi(y) < C minuj(y), 

\y\=r \y\=r 



maxfj(y) < C m\nvi{y) < C minify) < Cvi(0) < C 

\y\=r \y\=r \y\<r 



for < r < 1, and the claim follows. □ 
The next proposition is the analogue of proposition 4.3 of [32] and of proposition 1.4 of 
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Proposition 6.3. Let Xi — > x £ dM be an isolated blow-up point and assume that E4 — > 00 and 

n+2 
n-2 



€i — >• are given. Then pi — > ^| and, after passing to a subsequence 



Vi-U || C 2 ( o (0)) < e { . 



and 

R, 



log Mi * " 



Proof. Let R > and e > be given. Prom lemma IST21 we have < C. 

Therefore, by standard elliptic estimates there exists a subsequence of v,- t converging in Cf to a 
limit v which satisfies 

'Av + KvP = inlEr T , 
< = on <9M__ T if T < 00, 

u(0) = 1 and y = is a local maximum of u, 

Pi-l Pi- 1 

where T is the limit of a subsequence of Tj = M f 2 dist 9i (xj, dM) = M i 2 \xi\ = \y~i\. If T = 00 
then the proposition follows from the well known result of Caffarelli, Gidas, and Spruck (|14|). If 
T < 00 then the boundary converges to a hyperplane when i — > 00 by lemma I6.1} and the result 
follows from proposition lA.il □ 

The following lemma is analogous to lemma 2.1 of [23]. As in the the proof of [23] - - where 

pj-i 

they assume conformal flatness — the idea is to show that if the M i 2 dist ffi (xi , dM) does not stay 
bounded, then after rescaling the solutions we obtain an interior blow-up point, in which case the 
machinery of [25] can be applied (of course, in [23] they could not use [25j since such results had not 
yet been known, but they could still apply whatever was known about blow-up points in conformally 
flat manifolds without boundary; the idea here is similar). The proof does not require change to 
y-coordinates but we will keep track of the expression in y-coordinates for future use. 

Lemma 6.4. Let x% — > x S dM be an isolated simple blow-up point, with x^ € M. Then 

M~^ dist ffi (xi.dM) 

stays bounded. 

Proof. Let Xi be such that dist 9i (xj, dM) = dist Si (xj, Xj). The proof is by contradiction. Consider a 
subsequence such that 

M7*~ d\st gi {x h dM) = M**~\xi\ -> 00 
Pi- 1 

i.e., \y~i\ — > 00. Put Ti = M i 2 \xi\ = \yi\ and take normal coordinates at Xj. For \z\ < \xi\~ 1 o~ (where 
a comes from the definition of isolated blow-up point) define 

£i(z) = N^miN-** 1 z) 

where N- = 
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Pi- 1 

Notice that £j has the same form as V{ with iVj in place of Mj, so if (gi)(z)ki = (di)ki(^i 2 z ) we 
see that & satisfies 

[L-g£i + /T/f^f 4 = for \z\ < Ixil-V, 

1 Sft 6 = ^ & + ^ = on <9M, 

where /i(^) = 2 z). Since Xj is an isolated simple blow-up point for m we have U{(x) < 

2 2 2 

C|x| and then £j(z) < C\z\ "t- 1 . This, together with the fact that £(0) = {x^^- 1 Ui(0) = 



2 



2 

|xj| p i _1 Mj = — y oo as i — y oo implies that {0} is an interior isolated blow-up point for £j, 
hence we can use corollary 2.6 of [25] (with iVj instead of Mj and £, instead of Uj) and conclude that 
£i(0)£i -> w in C;oc(^-i ~~ where to > is the Euclidean Green's function for the Laplacian 

centered at (Euclidean because converges to the Euclidean metric) and = {z n > — 1}. It 
also follows that Bg^j = d v ~Xi + ^f^K^i = becomes in the limit J^- = on (9M" 1 . We have (see 
for instance [23]) 

w(z) = a\z\ 2 - n + A + 0(\z\), A>0. 

Define 

*M.«.V0 = ^£-5lV«l' + r(£)». 
Because is an interior blow-up point, we can use theorem 7.1 of [25] to get 

liminf / B(r,z,w,Vw) >0 

-^o y| 2 | =r 

and a direct computation gives 

jn— 1 1 



/" n — 2 

liminf / B(r, z,w,Vw) = A\S r ' 

r^o J\ z \ =r 2 



contradicting A > 0. □ 

Suppose Xi — )• x € <9M is an isolated simple blow-up point. In the notation of lemma [fT4l write 

Pi- 1 Pi- 1 Pi- 1 

Ti = M i 2 dist 9i (xj, <9M) = M i 2 |xj|. In y-coordinates this becomes Tj = 2 |xj| = 

By lemma 16.41 we cannot have T — > oo, and passing to a subsequence we have T j — > T < oo. 

Corresponding to the subsequence {IL} there is a subsequence {wj }. Applying proposition 16.31 to 

the {viA yields T = 0. Hence, we can hereafter assume that 

(6.2) l&HO. 

Proposition 6.5. Let Xj — > x £ <9M 6e an isolated simple blow-up point for the sequence {u{\ of 
positive solutions to h2.1\) . Then there exist constants C > 0, a > independent of i such that 

MiUi(x) > C _1 Gi(x,Xi), M**~ < |x| < a, 

MiUi(x) < C|x| 2 ~ n , |x| < a, 

where Gj(x, Xj) is f/ie Green's function for L gi centered at Xj iwii/i boundary condition B gi Gi(x, Xj) = 
on d'B a {xi). Moreover, after passing to a subsequence MjUj(x) — > G(x,x) in Cf oc (B a (x)\{x}) , where 
G(x, x) is the Green's function for L g centered at x with boundary condition B g G(x, x) = on 
d'B a (x). 
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Proof. The proof is an adaptation of the ideas from [32] and [23J using lemma IA.6I and proposition 
IO □ 

7. A FURTHER ESTIMATE ON dist ff< (Xi , 8M) . 

Let Xi — > x £ <9M be an isolated simple blow-up point. As in the boundaryless case, one of the 
main features of our proofs is the usual use of coordinates centered at the points Xj. If x; L £ M, lemma 
16.41 then gives an estimate for the distance of Xi to the boundary. Unfortunately this estimate is not 
enough for our purposes. In fact the results of sections 0] and [5] require the center of the coordinate 
system to be on the boundary. We therefore have to prove that we can pass to a subsequence such 
that Xi £ M. 

Proposition 7.1. Suppose Xi — > x is an isolated simple blow-up point. Then in boundary conformal 
normal coordinates at Xi, there exist constants a, C > 0, independent of i, such that 



U\(y) < Cei 



for every \y\ < oM i 



2 



Proof. The idea of the proof is as follows. Using the fact that the boundary is totally geodesic in 
boundary conformal normal coordinates, we can reflect all quantities across the boundary and then 
mimic the proofs of [32]. In order to simplify notation the index i will be dropped from all quantities 
when no confusion arises, and the metric g in y-coordinates will simply be denoted as g. Similarly / 
will be denoted by /. We will use Greek letters to denote indices running up to n and Latin letters 
for indices running up to n — 1, and write as usual y = (y', y n ). Let I = ere -1 . 

If yi € dM, take Fermi coordinates (z , . . . ,z n ) at yj. If € M then take Fermi coordinates 
{z\ . . . , z n ) at where yi G dM is the closest point to y^. Then in these coordinates g nn = 1 and 
g„,i = 0. Shrinking the domain if necessary, we can assume that the domain of definition of the Fermi 
coordinates contains the domain of definition of the boundary conformal normal coordinates. Define 
the extensions 

(7U -l , n, jg(z',Z n ), Z n >0 jv(z',Z n ), Z n >0 

(7.1) Q(z,z ) = < , and v(z,z ) = < 

K ' y ' \g(z',-z n ), z n <0 ( ' \v{z',-z n ), z n <0. 

Recall that in boundary conformal normal coordinates the mean curvature vanishes and hence the 
boundary condition for v is just a Neumann condition. Moreover the umbilicity of dM gives that 
the second fundamental form vanishes as well. Therefore the above extensions are C 2 , and are in 
fact smooth in the z' direction. Notice also that we are performing a change of coordinates to Fermi 
coordinates, but we are not making a conformal change of the metric, and hence the vanishing of k 
and Kij are still true in Fermi coordinates. Mimicking a standard one-dimensional argument then 
shows that d n (d^y) and d n {d^Jj) exist in the weak sense, so in particular the extensions are C 2,a . Of 
course, the extended metric also satisfies g nn = 1 and g ni = 0. 
A simple calculation shows that 

(7.2) R^z\z n ) = R g {z\-z n ), z n <0, 
and 

(7.3) Ajv(z', z n ) = A g v(z', -z n ), z n < 0. 
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Now extend the function / across the boundary by f(z', z n ) = f(z', —z n ) if z n < 0. Notice that / 
and Rg are C°> a . Combining (|7.2|) and (|7,3|) produces 

Lfv(z\ z n ) = L g v(z>, -z n ) = -Kf- S (z', -z n )v?{z', -z n ) = -KT\z', z n )v p {z', z n ) 

for z n < 0, i.e., the extended quantities also satisfy the equation. It follows that the extended 
equation holds in the original y-coordinates, 

(7.4) Lp}{y) + KT & v p {y) = in 5/(0), 

where 5/(0) is a full ball in R n , i.e., 5/(0) = {y G R n \ \y\ < I}. From det g = 1 + 0(r N ) in 5/(0) we 
obtain det g = 1 + 0{r N ) in Bi(0) as well. 

Now that the problem is defined in the full ball -B/(0), to prove the proposition, proceed with 
almost identical arguments as in the proofs of lemmas 5.1, 5.2 and 5.3 of [32J. There are, however, 
three differences that we now discuss. 

First, unlike in [32] the coefficients, of the PDE are not smooth. However, they are sufficiently 
regular to apply elliptic estimates. 

Second, we need the estimate v < CU. In [32j this arises from the fact that the blow-up is isolated 
simple. In the current situation it is not necessarily true that is an isolated simple blow-up point 
for v on Bi(0). Nevertheless, we will show that v(y) < CU(y) still holds for all y G Bi(0). Notice 
that we do not need to make an extension of U since it is a priori defined on the whole of W 1 . 

To see why this is the case, first notice that since yi is an isolated simple blow-up point for v on 
Bi(0), we have v < CU there. For p G 5/(0), let p G B/(0)\B/(0), be the reflected point. If ^ G dM 
then d-g(yi,p) = dj(yi,p), where c% means dist^. If yi ^ dM, then in y coordinates the boundary is 
given by a graph y n = F{y'), but F(y') — > as i — > (see (16.21) and lemma ISTTj) . which then implies 
dg{yup) = dg(yi,p) + o(l). Therefore 

v(p) = v(p) < CU(p) = C(l + dg( yi ,pf) 2 -^ < d(l + dg( yi ,p) 2 f-^ = dU(p), 

as desired. 

Finally, the third difference with [32J is that there, the scalar curvature satisfies R g = 0(r 2 ), which 
comes from the Taylor formula and properties of conformal normal coordinates. Here, since Rg is 
C°' a only, we avoid the Taylor expansion. Without Rg = 0(r 2 ) the proof in [32] yields a weaker 
estimate, but since we only need \v — U\(y) < Ce, the hypothesis Rg = 0(r 2 ) is not necessary. In 
[32] the better estimate \v — U\(y) < Ce s , with s > 1, is established. □ 

Remark 7.2. Observe that as in [32], the proof of proposition 17. II produces the estimate 8{ < Cei. 

In the proof of the next proposition, we retain the notation for the reflected quantities that appears 
in the proof of proposition 17.11 

Proposition 7.3. Under the same hypotheses of proposition p7T7| there exists a constant Co, inde- 
pendent of i, such that 

II v i - U \\ C 2, a (B h (o))— Co^i, 

where U = ere' 1 . 

Proof. It is sufficient to establish the desired estimate for uii =v% — U. We have 

L^Wi + biWi = Qi in B k (Q) 



26 DISCONZI AND KHURI 

with 

bi{y) =Kf 5i i (y), 
Vi u 

Qi{y) = {c{n)e 2 l R ai {e l y)U(y)+K(U^ - /~V>) + e? +1 0(\y\ N )\y\(l + \y\ 2 )- ^ 
Use ([7]) to find |&i(?/)| < c(l + |y|) -4 - Then the representation formula gives, for any and \y\ < ^, 

(7.5) Wi(y) = [ Gi(y,z)(biWi - Qi)(z)dz - f -J^-(y, z)wi(z)dS(z), 

JB k (o) JdB,.(p) w~gi 

where Gi is the Green's function for the conformal Laplacian with Dirichlet boundary condition. 
The proof is now similar to standard estimates for the Newtonian potential, and therefore we will 
only indicate the main steps (see for example [22J). 

First notice that unlike the Newtonian potential case, there is a boundary integral in the rep- 
resentation formula (|7.5p . Nevertheless, this boundary integral is easily estimated using standard 
properties of the Green's function and Vi < CU, since the singularities occur within the radius 4. 

For the interior integral, write ji = biWi — Qi. This quantity plays the role of the inhomogeneous 
term in potential theory. Therefore standard potential theoretic arguments yield 

C / \ 
( ? - 6 ) ^ 2 ^]a,%(0) < p( II 7i llco(IUo)) +™ ai5ii(0 )) 

where [-] a ,n is the Holder semi-norm on f2. Next, observe that by interpolation 
( 7 - 7 ) feto) ^ C (^L,B !i( o)K] Q ,5 ii(0 ) + ^4,^(0))" 

"2" "2" "2" "2" 

In order to estimate [wi] aB[ ^ the representation formula (17.51) may again be employed along with 

"2" 

standard properties of Gi and proposition 17.11 However, control of the boundary term relies on y 
staying away from the boundary, that is why we choose an estimate on Bu (0) (giving then a final 

~ n + 2 5. n + 2 

estimate on 5^(0)). Moreover, using remark 1731 and U™- 2 — f 'U Pi = U n ~ 2 0((\ log f\ + 1 log U\)5i), 

4 

it follows that 



(7-8) [Qi] a ,B h (o) < c[4[R 9 My)u{y)] a ^ (0) + [u~ - f •z/w] oJ 

+ £l [efO(\y\ N )\y\(l + W 2 -- 



2l a,B h (0) 



< Cei. 



Finally, the term || 7^ H^o^ (q)) * s estimated in a similar manner 

"2" 

II Ilc7°(5j. (0))- C^i- 



Combining this with ()7.6|) . fj7.7|) and (j7.8|) yields [D Wi} OL § l ( ) — C £ «- The remaining lower order 
terms of the C 2,a norm may be estimated in an analogous way. □ 
The analogous of the following result is already known for scalar- flat manifolds [2]. 
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Theorem 7.4. Suppose x% —> x is an isolated simple blow-up point. Then in boundary conformal 
normal coordinates at Xi, for all i sufficiently large and possibly after passing to a subsequence, we 
have Xi E dM. 

Proof. The proof is by contradiction. Therefore assume that Xi € dM occurs only for finitely many 
i. Hence passing to a subsequence, still denoted Xj, we can assume that 

(7.9) Xi € M for all i. 

Take boundary conformal normal coordinates at x% (see corollary 13. 10H . rescale all quantities to 
y coordinates as explained at the beginning of the text, and denote by € dM the closest point 
to yi, where yi is identified with the origin. The closure of the ball of radius \yi\ will be denoted by 
-Bigj(O). Furthermore, for any domain Q, denote by [-]i4- a ,n the C 1,a Holder semi-norm, and by 
the C 1 Holder semi-norm. 

Let Wi = Vi — U, then 

[ ' \iji-0\P S l^Ji+^, fii |(o)- 



As explained in section El the coordinates may be arranged such that q^- = g nn d n \g.. Observe 

that the boundary condition for Vi implies that d n Vi(yi) = 0, since the mean curvature vanishes. 
Notice also that we have Vvj(0) = VJ7(0) = 0. On the other hand a direct calculation gives 

(7.11) d n U{m) = (2 - n)(l + |yi| 2 r%. 
Hence (I7.10P becomes 

\m = \vi\ < ^^(i + l^l 2 ) f |y/Ni +A ^(o), 

since = (0, . . . , 0, yf). By (|6.2p . -^(l + \yi\ 2 )^ < C\ for a constant C\ independent of i, so 

(7.12) \yi\ < Cil&fiwi}^-^-^. 

By proposition 17.31 Wi converges to zero in C 2,a , so there exists a small r > 0, independent of i, such 
that the Taylor formula for Wi holds in B r (0) for all i. By (j6.2j) we can assume that B^(0) C B r (0). 
Therefore for any y 6 B\y.\(0), 

d k wi(y) = d k wi{0) + K^y 1 = Ki(y)y l , 

where we used V^j(O) = 0. The remainder term satisfies, for each I = 1, . . . , n, 

\Kl(y)\ < sup \V 2 w(z)\ <\\ Wi \\ C 2, a r B (0 ) )< C Si, 
*eB| S .|(o) 

where proposition 17.31 has been used. Hence \dkWi(y)\ < \TZi(y)y l \ < nCo£i\y\, and therefore 
(7-13) |(0 ) < nC Ei\yi\. 

Let O be a convex domain. The following inequality is standard (see e.g. 
(7.14) [«] 1+/ j, n < Ap a -P[u] 1+a , Q + Ap- p [u]i fl , 
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where the constant A depends only on the dimension, < /3 < a < 1, and p > is any positive num- 
ber. Also, using the mean value inequality, there is a constant A depending only on the dimension, 
such that \diu(jp) — diu{q)\ < A\p — q\[u]2,n = A\p — q\ a \p — q\ 1 ~ a [u]2,n- From this it follows that 

(7.15) [«]i+a,n < ^diam(fi) 1 - a [n] 2 , Q . 

Because the constants Co, C\,A and A do not depend on i, we can, with the help of (j6.2[) and the 
definition of E4 , choose i so large that 

(7-16) \m\ < ~ 

(7.17) £i < 1, 

<£ 1 

(7.18) max{C CiAA, nCoCxAje/ < -, 

(7.19) Code]-" < 1, 
where p > 1 is a large number chosen such that 

(7.20) ol———— < 1. 

V 

This is possible since a < 1; notice that p does not depend on i. 

Now fix an i = i (n, C , Ci,A, A) such that (f77T6l> - (j77T9]> hold. From (I7TT51) we have 

(7-21) Klu„CT<^K 



Jl oJi +Qj S| S . ,(0) - ^"^o^B,. ,(0)- 

Moreover the constants Co,C\,A and A do not depend on the choice of /3, as can be seen from 
the derivation of inequalities (j7. 12|) . ()7.14p . ()7.15j) . and the proof of proposition 17.11 Therefore the 
inequalities (f7TT2>(f?3T]) hold for any /3 such that < /3 < a. 

We are now in a position to prove the theorem. It will show by induction that 

(7-22) |& I<4 Q 

for all k = 0,1,2,3,... . Since £j < 1 this would imply \yi \ — so that Xj £ dM, contradicting 

For k = (|7.22p is true by (|7.16p . For k = 1, recall that d n Vi (yi ) = 0, and observe that 

\dnU(y io )\ = \d n (vi - U){y io )\ < Mi, B| -. ,(0) - 
Then (^TT]) . proposition O and (f77T9]) give 

\yi \ < Ci[wi ] 1B (0) < C Cie io = CoCie^e? < ef Q . 

' ,yi o' 

So assume that ()7.22p holds for some k > 1. Combining ()7.12|) and (|7.14p gives 

(7.23) \y l0 \ < d\y i0 f KJ^^-p ^ C iMVi f (p a ~^ [m } 1+a ^-^ + P^Kli,^- ^y) • 

Choose /3 = ^ (which is less than a by the choice of p). If we also choose p = e^ Q then (|7.23p 
becomes 

li«ol < CiA|^ h fc ^0) + £ 40 P K]i,^-^)J. 

By (|7.13p , (]7.2ip , proposition 17.11 the induction hypothesis (|7.22p and the fact that 

\WiJ„-B <ll Wi 



^J 2 ,B| Slo |(0) -II u '«0 IIC2,«(B |s . o| (0))' 
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we obtain 



, ^ — / fca+l-- -S+l+fca 

|£r*„| < max{C CiA^, nCodAje* U p + e p 



,2 



2max{C CiiVA nC CiiV}e> ° " 



< £ (fe+l)a £ l-«-p 



1-QJ-" 

where (f77T8|) has been employed. Finally, e io p < 1 by (fTTTj) and ([7720]) . □ 

8. Symmetry estimates 

In this section we derive sharp estimates for the behavior of solutions Uj in the neighborhood of 
an isolated simple blow-up point. The proofs are an adaptation of the results of |25] and we will 
often refer the reader to it for details. 

Throughout this section, let Xi —> x € dM be an isolated simple blow-up point. By theorem 
17.41 we can assume that x^ G dM. We will be using boundary conformal normal coordinates at Xi 
(see proposition 13. 5[) and rescale all the quantities to y-coordinates as explained at the beginning 
of the text. Notice that because in boundary conformal normal coordinates we have K g = in the 
neighborhood of the origin, the boundary condition becomes a Neumann condition. Moreover, since 
the boundary is umbilic we obtain that it is totally geodesic in the neighborhood of the origin. 

Also, by (|4.3p . proposition 15.11 gives that for \y\ < cre~ x we have, with zi = z E ., 



•1) 



Afi + n(n + 2)U— Zi = c(n) Yl=l didjH^U, for \y\ < ae^\ 
= J2?=i9 nl dik = e?0(\y'\ N (l + ly'l) 1 -™), on dM 



where (J2.9P has also been used. Notice that since d n U(y',0) = 0, in these coordinates U also satisfies 
the boundary condition 

Pjjj n - 1 

(8.2) — = Y,9 nl diU = efO{\y'\ N {l + \y'\) l ~ n ) on dM. 

v ' 9i i=i 

Proposition 8.1. Suppose Xi — > x is an isolated simple blow-up point. Then in boundary conformal 
normal coordinates at X{, there exist constants a, C > such that 

\vi-U-~ Zl \<C max {ef \H^\ 2 ( Xi ), e?~ 3 , Si} 

2<k<d-l 



for every \y\ < o~M 



vizi 

2 



Proof. Put Aj = maxi^i^. \vi — U — z%\ = \vi — U — Zi\(jji). Then as in the boundaryless case we get 
a stronger inequality if there exists a constant c such that \yi\ > cli for every i. In fact, using that 
x is an isolated simple blow-up point we get the inequality v < CU < C\y\ 2 ~ n , and using estimate 
(j2~9j) we get Aj = \vi - U - Zi\{yi) < C\yi\ 2 ~ n < e™~ 2 . Hence we can assume \yi\ < l - L 



2 ' 



If the proposition is false we have 
(8.3) 1 max {ef\H^\\ Xi )} 0, ^e^ 3 -► 0, 1* -+ 0. 

A; 2<k<d—l A; A ? - 



Define 
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Then |u>i(y)| < 1, and 
(8.4) 



' Lg.Wi + biWi = Qi in B h (0) 

Wi = 0(A -1 e^~ 2 ) on d+B k (0) 

& = A-^f 0(\y'\ N (l + W\) l - n ) on 5^(0), 



where (|8.ip . (|8.2|) and the boundary condition for Vi have been used; Qi and b{ are as in the bound- 
aryless case 

Vi-U - Zi 

1 n— 6 

6 — n n+2 ~ r 

+ 0(|^| 2 C/^) + K((C/ + 5;)^ " / + z) Pl ) 



0(\y\")\y\(l + 



|2\ — - 



and they satisfy the estimates (see 



\bi(y)\ <c(l + 



and 



|Q,(y)| <C-3-{ max {sflilWp^Xl + lyl) 2 ^ 2 -" 



+ £ r 3 (i + Mr 3 + m 4 (1+ ^o(|y|")| y |(i + \ y \ 

+ 5,(| log(E7 + 30| + | log fi\)(l + |y|)- n " 2 }. 

Let Gi be the Green's function for the conformal Laplacian with boundary condition d = on 
d + Bi i (0) and B gi Gi = = on d'Bi.(0). The representation formula then gives 



2\ — 2. 
Z ) 2 



Wi(y) 



r dG (v ri) 

Gi(y,rj)(biWi-Qi)(r])drj - / 10,(77) — 5— ^ — dS iv) 

Jd+ Bi . <r<\ OV5 



+ 



S'fl..(0) 



y, 



for 1 2/| < ^. The first two integrals are estimated as in the boundaryless case (see [25] ). For the third 
one we use (18.41) to find 



(8.5) 

Hence, 
(8.6) 



a's ; .(o) 



dv~ 9i 



< Ce 



n-2 



\m(y)\ < c[ (1 + M)- 2 + f max { £ f l^l 2 ^),^- 3 ,^} 

A,- 2<k<d— 1 
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It then follows from (18. 3D . (|8.6I) . and standard elliptic estimates that Wi is bounded in Cf oc and has 
a subsequence, still denoted Wi, converging to a limit Woo, which satisfies 

Awoo + n(n + 2)U^w oa = on 

= on R"- 1 

lim^^ Woc(y) = 0. 

Note that for the boundary condition we used that A~ 1 ef\y'\ N (1 + |y'|) 1_n < Ce™~ 3 for \y'\ < ae^ 1 . 
Lemma IA.5I then gives 

w = c (^^(y) + y ■ w) + ^ C A U - 

However u>j(0) = |Vu>j|(0) = implies u;oo(0) = |Vwoo|(0) = 0, from which we conclude that i^oo = 0. 
It then follows that \yi\ — > 00. This combined with (|8.3p contradicts (|8.6|) . as W{{yi) = 1. □ 

The proofs of the next two results are similar to those in |25j , making the necessary adaptations 
to the boundary case with ideas described in proposition 18.11 

Proposition 8.2. Under the hypotheses of proposition HOI 



5i<C max {ef \H^ \ 2 ( Xi ), e?~ 3 } 

~ 2<k<d-l l 1 ' ; 



Pi-l 

for every \y\ < oM i 2 



Proof. If the proposition is false we have 

(8.7) i max {ef l^p^l^O ier 3 

5; 2<k<d-l Si 1 



Hence from proposition 18.11 

\vi-U-Zi\(y) <C5i. 

Define 

1 , 

Wi = —{Vi - U - Zi), 

Oi 

and argue as in proposition 18.11 with 5{ replacing Aj, to obtain Wi — > Woo in Cf oc , where dnWoo = 
on IR™ -1 . Define ^(y) = 1! ^ L U{y) + y 3 djU(y). Now we argue as in [25], except possibly for the extra 
boundary terms 



T dwi r d$> 

w — — and / Wj 



d'B h (0) Jd'B k (0) dvgt 

But as before, §^ = efO{\y'\ N {\ + \y'\) l ~ n ), and a direct computation gives M- = efO{\y'\ N {l + 
|y'|) 1_n ), which is enough to handle the boundary integrals as in proposition 18.11 □ 

Proposition 8.3. Under the hypotheses of proposition HOI 

d-l 

|V m K -U- < CY,4 k \H (k) \ 2 (x t )(l + \y\)W~n-m + £ n-3 (1 + | y |)-l-»n 

k=2 

for every \y\ < ae^ 1 , m = 0, 1, 2. 
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_n— 3 



Proof. Arguing similarly to [25] with the necessary modifications as in propositions 18,11 and 18. 2} we 
obtain the result with m = 0. To obtain the result for the derivatives, we invoke standard elliptic 
theory, which gives the estimate provided that we can bound the C 1,a norm of d v ., {v{ — U — Zi) on 
the boundary. Since d v ~Vi = and d n Zi\ yn=0 = = d n U\ yn=0 , it is enough to show that 

n-l 

(8-8) || Y.^diizi + U) || o i,- (8 , fll . (0 ))< Csl 

i=i 

From (|2.9p . (|5. 1|) . properties of boundary conformal normal coordinates (in particular corollary [3 
and the explicit form of U we have 

n-l 

I Y,9 nl di{h + U)(y',0) < Cef\y'\ N {l + \y 



i\\l-n 



1=1 

which is bounded by Ce™~ 3 for \y'\ < cre^ 1 . 

Differentiating g nl di(zi + U) with respect to y k , k < n — 1, using again (|2,9p . (|5.ip . and 

properties of boundary conformal normal coordinates yields 

n-l 

\d k (Y,9 nl di(zi + U))(y',0)\ <Ce^ for \y'\ < os- 1 . 
i=i 

Differentiating again and repeating the argument gives (|8.8p . Now the pointwise estimate follows by 
standard arguments. □ 

9. Weyl vanishing 

In this section we will work mostly in x-coordinates and take boundary conformal normal coor- 
dinates at X{. In these coordinates, estimate (|2.8p and the estimate of proposition 18.31 become, for 
|x| < cr, 

„ n— 4 

n — 2 

(9-1) |V m z,(x)| < e t > \ h ^i + \x\)^ +2 ~ n ~ m 

\a\=4 jl 

n-2 ®Z3 n-2 

(9.2) \v^{ Ut - U£i - Zl ){ x )\<Ce^ ]T|tf {fc) | 2 (*0fe + M) 2fc+2 ~ n ~ m + V (e t + \x\)- m -\ 

k=2 

where both z\ and the sum with \H^\ 2 (xi) appear only when n > 6, and 

-l. 



z £i = Zi(x) = e 2 z i (e i x) 



n-2 



ti E; U = £T + a: 2 



Throughout this section it will be assumed that (M n ,g) is a Riemannian manifold of dimension 
3 < n < 24 with umbilic boundary. The index i will be dropped from all quantities in several 
estimates below. Note also that by theorem [73] we can assume that Xj € dM, therefore the boundary 
is given by dM = {x n = 0}. We will use the notation = {x G B p {xi) | x n > 0}, where p < a - 
of course, is the same as B p (0), but the first notation will be emphasized since it better suits the 
Pohozaev identity. Furthermore, the unit normal will be denoted by v = v g = v g . when no confusion 
arises, and vg will denote the Euclidean normal. 

We can now state one of the main estimates of the paper. 
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Proposition 9.1. Suppose 6 < n < 24 and that xi — > x £ dM is an isolated simple blow-up point. 
Then 

d n 

E N,a|V H |loge|^i <Ce n ~\ 

\a\=2i,j=l 

where 9^ = 1 if k = and 6^ = otherwise. 

Before giving a proof of proposition 19,11 some consequences are derived, in particular the Weyl 
vanishing theorem. 

Theorem 9.2. (Weyl vanishing) Let xi — > x be an isolated simple blow-up point and 6 < n < 24, 
then 

\V l g W g \\ Xl ) < Cer 6 - 2 '|log £j |-^, 

for every < I < f 1 ^] , where 9k = I if k = and Ok = otherwise. In particular \V g W g \ 2 (x) = 
for < I < . 

Proof. Proposition 19.11 gives the same estimate as in the boundaryless case, the argument then is 
similar (see [53]). □ 

Corollary 9.3. Under the same hypotheses of the Weyl vanishing theorem, 

iv>< -u- ziMi < cer 3 (i + ivir™- 1 

or, in x- coordinates 

n-2 



Proof. This is straightforward from proposition 18.31 and theorem 19.21 □ 

We now proceed with the proof of proposition 19.11 The proof will involve an application of the 
Pohozaev identity (|A.6P in a half ball . 

Write 4> = + x k dtu and 4> £ = ^^u e + x k dkU £ . In the proofs below extensive use will be 

made of the inequalities |V m u| < Ce^~\x\ n ~ m and |V m </>| < Ce^~\x\ ~ m , which follow from 
proposition 16.51 

Proof of proposition \9.1\ First it will be shown that there exists a constant C such that 

d n N n-1 

(9.3) c( E E N,^ 2H+1 + E E +z n ~ 



-2 

I ijiOi' I ^ ' ' ~~ I" ^ 

|a|=2i?=l |o'|=0i,J=l 



> - / 4> £ H in diZ £ dx + / c(n)4> £ u £ dijhijdx 
Jd'B+ Jb+ 

f n — 2 

+ / c(n)(4> £ z £ + u £ (— — z e + x k dkZ £ ))dijhijdx 
Jb+ 2 

f 11 

+ / c^^Usi-djiHijdiHu) + -djHijdiHu - -diHijdiHij)dx. 
Jb+ 1 4 

where by a' we mean derivatives along x' only. 
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Start with the Pohozaev identity (|A.6j) . On its left hand-side the integrals over the hemisphere 
"5+ 0°) are °f order e n ~ 2 and the boundary terms with x k v k vanish on d'B+ so the remaining term 
on the left hand side of ()A,6P is 

( - u + x k d k u)——dx'. 
&B+ 2 dug 

Since d v u = and g nn is bounded away from zero near the origin 



i=i 



therefore 



. n — 2 hn du , . 
(——u + x d k u)^—dx 
8B+ 2 dv s 



< C 



+ C 



.n-2 
( 1 

8B+ 2 

n - 2 

dB+ 



n-1 



+ x k d k u) 9 nl di(u -u £ - z £ )dx' 



1=1 

n-1 



II 



+ x k d k u) 9 nl di(u £ + z £ )dx' 



n-1 



Using (|4.2p . (|9.ip . (I9.2|) . and theorem 13. 4( we find that the above integrals are bounded by Ce 
and terms involving the umbilicity tensor. Now (|9.3p follows from (see [25]) 

- dijhij + dj{HijdiH*) ~ -djHijdiHu + ^H^H^ < C ^ E \ h ij,«\ 2 \x\ 2H + C\x\ n ~ 3 . 

\a\=2ij=l 

The next step is to show that, as in the standard case of a full ball, the first interior term on the 
right hand side of (|9.3p may be absorbed into the error. To see this observe that theorem 13.41 implies 

n-1 

,U £ (J2 9 H nO + d nHnn) + 0(e n ~ 2 ) 
3=1 

N n-1 

= °(E E^- 

|a'|=0i,i=l 

Therefore after an integration by parts (|9.3|) becomes 

d n N n-1 



Bl 



d'B2 



\a'\+l +£ n-2 



(9.4) 



\ a \=2ij=l 



2 £ 2|a|+l 



+ E E i r w 



As 



W\+l +e n-2 



> 



\a'\=Oi,j=l 

(c(n)4> £ u £ H in diHii - (j) £ H in diZ £ ) 



d'B2 



-2 fc(n)u £ z £ (l + ^x k d k )dijhij + [ c(n)4> £ u £ (^djHijdiHu - ^HijOiHij) . 

J Bp - J Bp 

Furthermore the boundary integral on the right hand side of (|9.4p may be absorbed into the error 
term with the help of theorem 13.4^ 

. N n-1 

/ (cWfeueHtodiHa - 4> £ H in diZ £ ) = 0( V V \T lha ,\e\ a '\ +1 + e n ~ 2 



\a'\=0i,j=l 
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The remaining interior integrals are the same as those that appear in the original Weyl vanishing 
proof [25J, except that the domain of integration is a half ball instead of the full ball. At this point we 
may follow the original proof to obtain the desired conclusion, as long as the following two facts hold: 
(i) the necessary integration by parts may be performed with the extra boundary integrals (along 
d'Bp) being absorbed into the error, (ii) an orthogonality condition among harmonic polynomials 
holds on the half ball. 

An inspection of the original proof shows that (i) is valid, since any integrand along d'B^ will 
contain quantities that either appear in theorem 13.41 (and thus may be estimated by the umbilicity 
tensor) or involve d n z £ — which vanishes by proposition l5.11 Furthermore consider the decomposition 
(|5.ip . then in the notation of |25j 

(H q )\f =VTomdiPk-2 q \x\ 2q+2 )- 
In section[5]it was shown that d n pk-2q\Qi vanishes, therefore it follows that (H q )i n , d n (H q ) nn , and 

d n (H q )ij can be estimated in terms of the umbilicity tensor. This implies that the corresponding 
elements of Wij can also be estimated in terms of the umbilicity tensor (since the corresponding 
elements of Hij have this property by consequence of theorem 13. 4p . Hence ((H q )ij, Wij) can be 
absorbed into the error, where the inner product is taken over the half sphere. Similarly 

/ (I- k)pip k = / (PkdnPi ~ PidnPk) = 0, 

JS^ip) Jd'B+ 

so that pi _L pk, I 7^ k, that is, (ii) is valid. This finishes the proof of proposition 19.11 

10. Sign restriction 

Define 

P'(r, w) = / {^— W 1T- + xkd k w a -x k v k 8 \Vw\ 2 )ds. 

JdB+(x t ) 2 dv s dv s 2 

Proposition 10.1. (Sign restriction) Let Xj —> x be an isolated simple blow-up point and assume 
that 3 < n < 24. If MiUi{x) —> w away from the origin then 

liminf P'(r, w) > 0. 
r-»0 



Proof. Define 

P(r, Ui ) = I {^Uip- + Jdnutp- - W|Vtn| a + -^—K(x)x k uluf +1 )ds. 

If r is sufficiently small, the Pohozaev identity (proposition IA.7D gives 

P(r, Ui) > - I (c-r^Ui + x k d k Ui) ((g 1 / - 5 h )dijUi + dig\ 3 djU^dx 

JB+{xi) 1 

C o u i + x k d k Ui)R gi Uidx. 
B?{xi) 2 



Notice that 

' L>} 1 r,) 



(10.1) / x k d k UiR gi Uidx 

JB+(xi) 
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= -/ {x k d k R 9i +nR 9i )u 2 dx - x k d k u i Rg i u i dx + x k v k h R 9i u 2 . 

JB+(xi) JB+(x l ) J8B+(xi) 

Since x k v k = on d' B^~(xi) and v k = x k jr on d + B+(xi) we obtain 

/ x k d k UiR gi Uidx = -- I {x k d k R 9i +nR 9i )u 2 dx + T - [ R 9t u 2 ds, 

JB+(xi) 1 JB+(x,) 1 Jd+B+(xi) 

SO 

r n — 2 

(10.2) c(n) / ( m + x k d k Ui)R gi Uidx 

JB+{ Xi ) 1 

= -c(n)/ (-x k d k R 9i + Rg^ujdx + c{n) r - j R 9l u 2 ds, 

JB+(x t ) z 1 Jd+B+(xi) 



and then 



P(r,Ui) > - f (^-^Ui + x k d k u i )((g l j J - 5 h )dijUi + dig l /djUi)dx 

JB+(x t ) 1 

~ c(n) / {\x k d k R 9i +R 9i )u 2 dx + c{n) r - R m u 2 ds 
JB+(xi) 2 2 Jd+B+ixi) 

= Ai(r) + c{n) r - I R 9i u 2 ds 

z Jd+B+{xi) 

where Ai(r) is defined by the above equality. Now observe that fgB+( Xi ) K{x)M 2 x k v k u p i %+l — > 0. In 
fact, the integral over d'B+(xi) vanishes as x k v k = there. On d + B^~(xi) we have x k u k = r, hence, 
using the equation satisfied by Uj produces 

/ K{x)M 2 uf +1 = - [ MiuLgMiUi -> - f wL g w = 0. 

Jd+B+{xi) Jd+B+{xi) Jd+B+{xi) 

Therefore M 2 P(r,Ui) — > P'(r,w), so 



P'(r,w) = lim M?P(r,Ui) > lim M 2 Ai(r) + lim c(n)- / R g .(MiUi) 2 ds 

i-Kso i->co i^co 2 Ja+B+(xi) 

= lim Mf Ai(r) + c(n)~ / R 9l w 2 ds. 
t ^ 00 2 Ja+B+ixi) 



ld+B+( Xi ) 

We now proceed to analyze M 2 Ai(r), noticing that since theorem 19.21 and corollary 19.31 give the 
same estimates as in the boundaryless case, the same analysis can be carried out, except for an extra 
boundary term that appears in Ai(r) when integration by parts is performed, where 



"(x<) 

{^-(uet + x k d k (u Ei + z Ei ))(A 9i - A s )(u £t + z £i )dx. 

i) 2 



(10.3) Mr) = -c(n) [ {\x k d k R 9l + R 9i )(u £i + z £ fdx 

JB+(x,) * 

-f 

Corollary [93] implies that E 2 ~ n \Ai{r) - i;(r)| < Cr, so lim^oo e 2 i ~ n {A i {r) - Ai(r)) > -Cr. Notice 

2 

that since M{ = e i Pl ~ 1 and 2 - n we can replace M 2 by s 2 ~ n and obtain 

(10.4) P'(r,w) > -Cr + c(n)- [ R 9i w 2 ds + lim e 2 ~ n Ai(r). 

2 Ja+B+ixi) 



COMPACTNESS AND NON-COMPACTNESS FOR THE YAMABE PROBLEM 



37 



Using the symmetries of u £i 



(10.5) 



,n-2 



(-s-K + + x K d k (u £i + z ei )){A gi - A 5 )(u £i + z £i )dx 



,n-2 



+ 



r 

dB+(xi) 



(^—z £i + x d k z £i )(A 9i - A s )z £i dx 



,n — 2 u dz £j dz £ . . , 

(— u £i+X k d k u £i )(^--^)ds 

z JLu £i _dp ±)d ^ 



du, 



where L 



n-2 



,, + x k d k . The integrals over d + B^~(xi) vanish by properties of normal coordinates, so 

dz 



consider the integrals over d'B+(xi). Observe that = = d n z £i by proposition 15.11 Then using 
(|4.2p . the definition of u £i , and (|9.ip . we obtain 



(10.6) 



,ra-2 



-u £i + x k d k u £ 



d'B+{xi) 



dz £ 
) — - 
dv 



-ds 



,n-2 



n-l 



-U E 



d'B+(xi) 



< 

l8'B+{ Xi ) 

< CeT 2 r. 



£j + X 



+ x k d k u £i ) s ^gfd l z £t ds 
i=i 

- n {ei + \x'\f- n - l \x'\ N dx' 



For the other boundary integral notice that 
dLu r . dLu F , 



n-l 



dv gi dv s 



{-$°d a + d n )Lu £i = (-gr + l)Lu £i - Y,9i%Lu £i . 



1=1 



Since |VLu £i | < e i 2 (e$ + |x|) n , using (|9.ip . theorem l4.lt theorem |3.4| and (|4.2p . it follows that 

dLu £ dLu F 



(10.7) 



9^ 



< Cer 2 r. 



Combining (flTX3D . ([Tail . (flQ3|) . (fT03|) . and (flOTD yields 



(10.8) 



P'(r,w) > -Cr 



+ c(n) 



,n- 2 



z £i + a; d k z £i )(A gt - A s )z £i dx 



R„,w ds 



d+B?( Xl ) 



c(n) lim e 2 "™ / d* fc d fcJ R 9i + i^)(« £j 



We can now proceed as in the boundaryless case. The first integral on the right hand side of (110.81) 
as well as § Jd+B + ( x ) Rgi w2 ds are estimated using theorem 19.21 Theorem 19.21 and corollary 19.31 may 
be used to estimate f B + t yX .^{\x k d k R gi + R gi )z^.dx. Finally the estimate of proposition 19.11 is used to 
handle ! B +i x .){\x k d k R 9i + R 9i )(u 2 e . + 2u £t z £i )dx. □ 
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11. Blow-up set 

In this section we show that the set of blow-up points is finite and consists only of isolated 
simple blow-up points. The proofs are very similar to the boundary less case (|25j) and the locally 
conformally flat case with boundary (|23j). and therefore we will go through them rather quickly, 
indicating the necessary modifications. 

The following proposition is proven in |23j (proposition 1.1, see also |25} 132]). 

Proposition 11.1. Given 5 > sufficiently small and R > sufficiently large, there exists a 
constant C = C(5, R) > such that if u is a positive solution of \2. 1\) with maxu > C, then there 
exists {xi, . . . ,xj\f} C M, N = N(u) > 1, where ^| — p < 6 and each xi is a local maximum of u 
such that: 



1) {B Ti (xi)}f =l is a disjoint collection if ri = Ru(xi 



p-i 

2 



2) in normal coordinates centered at Xi 

|| Uiixi^uiuiixi)-^ y) - U{y) \\cp(B R (o))< $ 

where y = n(xj) £ 2~x, 

2 

3) u(x) < Cd g (x, {xi, . . . , x n }) p- 1 for all x € M and 

2 , 

d g (xi,Xj)p~ 1 u(xj) > C 

for Xi / xj. 

Lemma 11.2. Let Xi x be an isolated blow-up point for the sequence {ui} of positive solutions of 
\2.1\) . Then x is an isolated simple blow-up point. 



Proof. We argue as in [25j to obtain a subsequence Wi such that 

Wi (0) Wi (y) -+ h{y) = a\y\ 2 - n + b(y) in Q 2 oc (M^\{0}), 

where b(y) is harmonic in and satisfies d n b = on IR™^ 1 . Therefore, extending b across R n_1 and 
using Liouville's theorem shows that a = b > 0. Arguing as in [25J this leads to a contradiction with 
proposition 110.11 □ 

Proposition 11.3. Let 5, R, u, C(5,R), and {xi, . . . ,xn} be as in proposition \11.1[ If 6 is 

sufficiently small and R sufficiently large, then there exists a constant C(5,R) > such that if 
maxjv/ u > C then d g (xj,xi) > C for all 1 < j ^ I < N . 

Proof. Again we argue as in [25J, making the necessary modifications along the lines of [23J as in 
lemma rjOJ □ 

The following is an immediate consequence. 

Corollary 11.4. Let {ui} be a sequence of solutions of \2. 1\) with rn.ax.MUi — > oo. Then pi — > ^| 
and the set of blow-up points is finite and consists only of isolated simple blow-up points. 

12. Compactness 

Now that we have the Weyl vanishing theorem and sign restriction, the remaining arguments for 
the proof of theorems 11.11 and 11.21 are similar to those of the boundaryless case. In fact, the results 
of this section will be an adaptation of [25l[8l[T7], and therefore as in section [Til we will go through 
the proofs very briefly. 
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Proof of theorem From the results of section [TT] pi — > and there exists a finite number 

N > of isolated simple blow-up points x ^ y x^ ^j...j3/- y x { K If none of the X£ belong to 
the boundary then the compactness result follows from [25], so assume that at least one of them 
belongs to dM. It may also be assumed without loss of generality that X£ € dM, £ = 1, . . . , N — k 
and X£ $l dM, £ = N — k + 1, . . . , N, for some k < N — 1. Furthermore let 

Ui (xf ) ) = mm{ui ) , . . . , ui {xf'^ ) } 

for all i. 

Set Wi = Ui(x^)ui. A standard estimate gives that away from the blow-up points Wi — > 

^2f=i a jG x u) j where dj > 0, a\ > and is the Green's function for the conformal Lapla- 

cian with singularity at x^> . Now argue as in [H] (see [H [T7] as well) to obtain the asymptotic 
expansion 

n-2 

(12.1) G(x,x^) = \x\ 2 - n (l+ fa) +A + 0(|x|log|x|), 

fc=d+l 

where G = G s (i), ipk are homogeneous polynomials of degree k and ^4 is a constant. The sum 
between parenthesis starts at k = d + 1 because hij >a (x) = at a blow-up point x € 9M, by the 
Weyl vanishing theorem. We remark that when the boundary is not umbilic an extra singular term 
appears in this expansion (see e.g. |26j). Also notice that standard properties of conformal normal 
coordinates, theorem 13.41 and the umbilicity of the boundary, imply that fgn-i dijHij = 0, from 

which it follows that 



(12.2) 




Now put g = G n ~ 2 g. Then (M\{x^},g) is scalar flat and its boundary is totally geodesic. If we 
introduce the asymptotic coordinates y = |x| _2 x, then the expansion ()12.ip and the Weyl vanishing 
theorem give %j = 5ij + 0(\y\~ d ~ l ). Therefore the doubling of (M\{x^ } ,g) is asymptotically flat 
and has a well defined ADM mass (|31j. compare also with [8]). 

The rest of the argument now is standard. The positive mass theorem (see remark below) along 
with (|12.2p and the Weyl vanishing give that A > (as in [25], using the hypothesis that the 
manifold is not conformally equivalent to the round hemisphere we can rule out the A = case). 
This contradicts the sign restriction of theorem HU.li finishing the proof. 

Remark 12.1. 

1) Strictly speaking, we did not show how to prove a positive mass theorem (PMT) for manifolds 
with boundary, as the mass of such manifolds was never defined. What is referred to as the PMT 
for manifolds with boundary is actually the statement that the constant term in the asymptotic 
expansion of the Green's function is non-negative, which in turn is implied by the positivity of the 
mass of the doubled manifold (see [32]). 

2) The PMT is known to hold up to dimension 7 [4U1 W2\ and in arbitrary dimensions if the 
manifold is spin [45 \ I27| . Therefore, our result for n > 8 in the case of non-spin manifolds is true 
provided that the PMT holds under such hypotheses. 



Proof of theorem This follows from lemma 111.21 and theorem 19.21 
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13. Blow-up of solutions for n > 25 

In this section we prove theorem 11.41 We assume n > 25 throughout. As we mention in the 
introduction, the proof relies heavily on the constructions of Brendle [7J and Brendle and Marques 
[9], and we refer the reader to them on several occasions. 

We start collecting facts from [7J 19] that will be of direct use in our proof. Their main results is 

Theorem 13.1. (Brendle and Marques, ]9\) Assume that n > 25. Then there exists a metric g 
on S n (of class C°°) and a sequence of positive functions m £ C°°(S n ) with the following properties: 

(a) g is a small perturbation of the round metric go which is not conformally flat, and g = go 
near and beyond the equator, 

(b) for each i, U{ is a solution of the Yamabe equation 

n+2 

Lgm + Ku"" 2 = 0, 

where K = n(n — 2) is a positive constant, 

(c) E g (ui) < Y(S n ) for all i G N, and E g {ui) — > Y(S n ) as i — > oo, where E g (ui) is the Yamabe 
energy of Ui and Y(S n ) is the Yamabe invariant of the round sphere, 

(d) sup 5 n Ui — ?> oo as i — > oo. 

The scalar curvature of the metric g satisfies 

(13.1) R g > c> 0. 

for some constant c, since g is a small perturbation of the round metric. In particular this guarantees 
the coercivity of L g , which allows us to use the the C°-blow-up theory developed by Druet, Hebey 
and Robert [16]. From their results and estimate (c) of theorem 113.11 it then follows (theorem 5.2 
of [16], see also discussion at the end of section 5.1) that Ui has only one blow-up point, and it is 
apparent from [7J[9] that this is the south pole (from the point of view of stereographic projection). 
Moreover, up to a subsequence the following estimate holds (again theorem 5.2 of |16| ) 

(13.2) Q^Ue^xiix) < Ui{x) < Qu ei)Xl (x), 

2 2 

for some constant Q > 1 independent of i and for all x 6 S n ; here Ei = (sup SlI Uj) ™- 2 = Ui{xi) n -' 2 , 

n— 2 2-n 

and u £i)Xi = e i 2 (e? + \x — Xi\ 2 )~ , \x - Xi\ = dist 9 (x,Xj). 

Consider now the south hemisphere , which we identify with the unit ball in R n via stereographic 
projection. Since g = go on a neighborhood <9£™, we have that 95™ is totally geodesic, and in 
particular B g = d Vg . Combining f)13.2|) with the Harnack inequality implies that away from the south 

2-n 

pole, e- 2 Ui converges in C 2 to a positive Green's function for the conformal Laplacian (possibly 
after passing to a subsequence). We claim that for large i 

dui 

13.3 ^ < 0. 

4 

To see this, denote by 5 the Euclidean metric so that go = AU n ~ 2 5. Let G go and G$ be the corre- 

n — 2 -, 

sponding Green's functions with singularity at zero. Their relation is given by G go =4 2 U Gg. 
Using (|A.3|) and the fact that the mean curvature of dS™ vanishes , we have 

^=B go G ao =4r*?ir£*B s G s <0 
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on dS™, where the inequality follows by direct calculation. Therefore -^f- < by theorem 113. lf a). 
so that pA3|) holds. 
We conclude that 

I L gUi + Kur 2 = , in 52, 
[B gUi < 0, on dSl. 



That is, Ui is a sub-solution of the boundary value problem 
(13.4) 



n + 2 

L g u + iTu^=0, in 5™, 
B 9 u = 0, on 052. 



Actual solutions to (|13.4|) will be constructed by finding appropriate super-solutions. The super- 
solutions will satisfy the equation with a different constant K, and this will require a slight modifi- 
cation of the standard sub-super-solutions argument. 

Theorem 13.2. For all sufficiently large i there exists a solution V\ of {13.$ satisfying < Uj. In 
particular sup S n Vi — > oo as i — > oo. 

Proof. Because of (|13.ip . we can choose 5 > so small that 

n+2 n + 2 

L g 5 + K6*=5 = -c{n)R g S + Kb~ < 0. 
Put Wi = AiS, where Ai > 1 is a constant chosen so large that 

(13.5) Ui < A^, 
and 

(13.6) -u;- 2 - < 
By the choice of 5 

(13.7) 

where -fQ = n ~ 2 K. So u?j is a super-solution of the problem with constant Ki. As pointed out 
before, (L g ,B g ) is invertible and therefore the operators T and Pi given by 




(13.8) Tz = t^ 
and 

(13.9) P iW = Pi & 



n+2 

L g t = -Kz^, in 5™, 
B g t = 0, ondS'l. 



~ n+2 

tfiii;^, in 5™, 
on 952, 




are well defined. By the maximum principle T and P\ are monotone in the sense that z$ < z\ => 
7"zo < Tzi, and analogously for Pj. 

Now we put u\ = Ui, iv® = uii and define inductively u^ +1 = Tu\ and wl +l = P%w\. Since Ui is 
a sub-solution we obtain u\ < u\ and inductively u\ < u^ +1 . Analogously w\ > since is a 

super-solution. 
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We have u® < by (113.5H . and claim that u\ < w;^ for every t (the difference from the standard 
sub-super-solutions argument is that the equations involved in the definition of T and Pi are not 
exactly the same due to the different constants K and Ki). The difference satisfies 



(13.10) 



L 9 (4 +1 - = -(K(u\)^ - Kiiwj)^), in S r l, 

B g {u<: +1 -w^ 1 ) = 0, on dSl. 



In order to apply the maximum principle we need the right hand side of (113, 10p to be non-negative. 
To show this, recall the definition of Wi and Ki, use the monotonicity of the sequences u\ and wf, as 
well as (|13.6p to find 

- n+2 ~ „ n + 2 „ n + 2 ~ „ n + 2 2±2 n + 2 

K(uj)^=2 - Ki{wl)^ > K(u°)^=i - Ki(w°)^ = Ku?~ 2 - KArf^ < 0, 



It follows that u\ < wf. 

Now a standard argument produces the desired solution uf of (|13.4[) such that u% < uf . The 
proof also yields a wf 3 solving (|13.4p with Ki in place of K, and such that and wf° < W{, but this is 
not the solution we are looking for due to the different z-dependent constant Ki. □ 

14. Leray-Schauder degree of solutions 

Here we discuss some consequences of theorem 1 1.11 Throughout this section we assume 3 < n < 24. 
The results here are very similar to the cases of manifolds without boundary and locally conformally 
flat with boundary, so we refer the reader to |25j and [23] for details. 

As we pointed out in the introduction, one obvious consequence of theorem 11.11 is to give an 
alternative proof of the solution to the Yamabe problem. This follows from the fact that standard 
variational methods can be used to give solutions to the subcritical problem 



(14.1) 



Lgii + KvP = 0, in M, 
B g u = 0, on dM, 



with 1 < p < More generally, the compactness theorem allows us to compute the total Leray- 
Schauder degree of all solutions to equation (|14.ip . and to obtain more refined existence theorems 
which we now discuss. 

Without loss of generality we can assume that R g > and n g = 0. Then we can write (|14.1j) as 



(14.2) 
where 



L g u + E(u)u p = 0, in M, 
|| = 0, ondM, 



E{u)= / (| V g u\ 2 + c{n)R g u 2 )dV g 

J M 



is the energy of u (there is no boundary term since K g = 0). Notice that the Neumann problem for 
the conformal Laplacian is invertible in that R g > 0. Defining 

Q A = {u£ C 2 ' a (M) | || u || C 2,a (M) < A, u > A' 1 } 

we obtain a map F p : — > C 2,a (M) given by F p (u) = u + L~ 1 (E(u)u p ). 

From elliptic theory, we know that the map u h- > L~ 1 (E(u)u p ) is a compact map from Q\ into 
C 2,a (M). Thus F p is of the form 7+compact, and we may define the Leray-Schauder degree (see 
of F p in the region f^A with respect to G C 2,a (M), denoted by deg(F p , Q\, 0), provided that 
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^ Fp(dCl\). The degree is an integer which counts with multiplicity the number of times that the 
value is taken on by the map F p . Notice that F p (u) = if and only if u is a solution of (|14.2p . 
Furthermore, the homotopy invariance of the degree tells us that deg(F p , Q\, 0) is constant for all 
P e [1' n~^§] P rov ided that ^ F p (8Qa) for allp G [1, Sirf]- Moreover, in the linear case when p = 1, it 
is not difficult to calculate, by an argument similar to what is done in [38J, that deg(7i, Qa, 0) = —1 
for all A sufficiently large. Therefore, theorem 11.11 allows us to calculate the degree for all p € [1, !rz§]- 
Since it follows from the a priori estimates we derived that does not belong to F p (dQ\), we obtain 



Theorem 14.1. Let (M n ,g) satisfy the assumptions of theorem \l.l\ Then for all A sufficiently large 
2^§], we have deg(F p , 



and all p € [1, ^jf], we have deg(F p , f^A, 0) = —1. 



In the case that all solutions of the Yamabe problem are nondegenerate, our previous results assert 
that there will be a finite number of solutions of the variational problem. Moreover, the strong Morse 
inequalities will hold for the Yamabe problem since these inequalities hold for subcritical equations, 
and theorem 11.11 shows that all critical points converge as p — >• ^| . It follows that 

A 

(-1) A < £(-l) A -"C„, A = 0,1,2,... 

/Lt=0 

where denotes the number of solutions of Morse index /i. Since there is a finite number of 
solutions, we then obtain: 

Theorem 14.2. Let (M n ,g) satisfy the assumptions of theorem \ and suppose that all critical 
points in [g] are nondegenerate. Then there is a finite number of critical points gx> - - - >9k> an d we 
have 

i = £(-i) /te) 

where I(gj) denotes the Morse index of the variational problem with volume constraint. 

Appendix A. Auxiliary results 

In this section we state several auxiliary results that are either well known or slight modifications 
of standard results. Therefore proofs, when provided, will be rather short. 

The following proposition is analogous to a well known theorem of Caffarelli, Gidas, and Spruck 
(CD): 

Proposition A.l. Let T > and M™ T = {y € R n | y n > —T}. Consider the problem 



Au + Ku p = 0, u > in 

& = ™ 

u(Q) = 1, is a local maximum of u, 
where p £ (1, ^r§]. If p < then this problem has no solution. If p = ^| then 

n-2 
2 



u(x',x n ) = {——— ™) = 17(a) 

VI + \{x',x n )\ z J 

in which case T = necessarily. 

Proof. |28] (see also the proof of proposition 2.4 in |20| . and |23| p. 498). □ 
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Now we recall some transformation laws. 
Proposition A. 2. Let (M,g) be a Riemannian manifold with boundary and<p > a smooth function. 

4 

Let g = 4> n - 2 g, then 

1 ri+2 

(A.l) Lg{(t>- l U) = <p-—lLgU 



n+2 



(A.2) R g = -c(n)-^ —L g<1 

(A.3) B- g {(j)- l u) = (p"^B g u 



4-n 



(A.4) kij = 0n-a Kij + ^7,<P n - 2 g^9ij 

2 n 

(A.5) k = -^^4> n - 2 B g <j> 

where quantities with ~ refer to the metric g, and k are the second fundamental form and the 
mean curvature, respectively. 

Proof. Direct calculation (see |25[ [T71 [32] for example) . □ 

Proposition A.3. Up to a conformal change we can assume that in small balls the scalar curvature 
is positive and that the mean curvature of dM vanishes. 

Proof. The idea of the proof is to perform two conformal changes on the metric, one to produce a 
metric with zero mean curvature and a further one to achieve positive scalar curvature. Denote by 
4>i > 0, the first eigenfunction of the conformal Laplacian with boundary condition B g (p\ = 0, i.e., 

(Lgfa + Ai0i = 0, in M, 
[Bgfa = 0, on dM. 

See [T7] for the existence of c/>i; the fact that (pi > follows from a standard calculus of variation 

4 

argument. By transformation law (|A.5|) . the metric g\ = (fti~ 2 g has zero mean curvature. 

Now let xq E dM and consider a small ball B2s{xq) near the boundary. Denote by tpi > 0, the 
first eigenfunction of the Laplacian A gi with the boundary condition as below: 

' A 9l *0i + mip! = 0, in B 2 s(xo), 
ipi = 0, on d + B 2 s(x ), 



B gi ^ = ^ = 0, ond'B 2S (x ). 



The existence and positivity of ip\ again follows from a standard calculus of variations argument. 

4 

Consider the metric g = ipi~ 2 gi on B2s{xq). Then from (1A.2I) . 

R-g = -c{n)- x ^—Lg^ x = -c(n)"Vr 2 (A^i - R gi ipi). 
Since fii — > oo as 5 — > we can choose 5 > so small that 

A gi ipi - Rg^! = -H\il>i - R gi ipi < 0, 
and therefore Rg > on Bs(xo). Notice that shrinking B 2 s(xo) does not affect R gi as <p\ is defined 

O _ n 

on the whole of M. Finally, the mean curvature for g is k = -^zz^ ™- 2 B gi vp = 0. □ 
The next result immediately follows. 
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Corollary A. 4. Up to a conformal change the maximum principle holds for the conformal Laplacian 
in small balls. More precisely, if L g u > in B a {xo), u > 0, then there exists a constant C > ; 
independent of u, such that sup B(r ^ X0 ^ u < C supg Bcr i X0 \ u, provided a is small enough. 

Lemma A. 5. Let ip be a solution of 



Ai/j + n(n + 2)U"- 2 ip = 0, in 



d± - n 
Qyn — U, 



on 



pn— 1 



Then it takes the form 



n-l 



Mv) = co (^f 7 + y ■ vu) + £ Cj djU, 



3=1 



for some constants cq, . . . , c n -\. 



Proof. Since d n ip = on M n , we can make a C reflection across R n and then the result follows 

□ 



from [15] . 

The following is a Harnack-type inequality. 

Lemma A. 6. Let Xi x be an isolated blow-up point and assume that f is sufficiently small. Then 
for all r such that < r < f we have 

sup Ui < C inf itj, 
B r {xi)\B r/2 (xi) B r (pH)\B r/2 {xi) 

for some constant C independent of i and r. 

Proof. It follows from a combination of lemma A.l of [23], the Harnack inequality, and the definition 
of isolated blow-up points. □ 

Proposition A. 7. (Pohozaev identity) Let u > be a solution of L g u + K ff S u p = on = {x € 
B p (0) | x n > 0}. Then 



(A.6) 



dB2 



^-l u + x k d k u)^- - \x k 4\V Q u\ 2 + J—K(x)x k iy k u p+1 )da 



p + l 



+ x k d k u) ((g* - 5^)d ijU + d jg ^d iU )dx 

B~X * 



f n — 2 1 

+ / c(n)( u + x k d k u)Rudx H 

Jb+ 2 p+l 



x k d k K{x)u p+l d. 



x 



B1, 



+ 



n 



n-2 



K(x)u p+1 dx, 



,p+l 2 , JB} 
where quantities with q refer to the Euclidean metric and K{x) = Kf~ s (x). 
Proof. Standard integration by parts argument. 



□ 
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